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THE BOULDER BATHYLITH OF MONTANA* 
ADOLPH KNOPF 


ABSTRACT. The Boulder bathylith is a large plutonic mass extending 70 miles along 
the Continental Divide of the Rocky Mountains from Mullan Pass on the north to 20 
miles south of Butte, Montana. The present investigation was begun in 1939 but was 
interrupted by the war years, and has concerned itself mainly with the northern portion 
of the bathylith and its enveloping rocks. 

The rocks invaded by the bathylith range from Beltian (Precambrian) age to late 
Cretaceous. The youngest prebathylithic rocks consist of andesitic and potassium-rock 
silicie lavas (dellenites and toscanites) and associated pyroclastics, In part the volcanics 
were submarine of Niobrara age and in part were erupted subaerially to form welded 
tuffs, probably in late Montanan time. Shortly after the final volcanism, the first severe 
Laramide orogenic crisis in this part of Montana set in. 

On the basis of a recent potassium-argon age determination, which gave 87 million 
years as the most probable age of the granodiorite of the Boulder bathylith, it is con- 
cluded that the bathylith was emplaced late in Cretaceous time. 

The Boulder bathylith has hitherto been considered to be a one-magma intrusion, 
but like other large plutonic masses it proves to be of composite construction, The order 
of intrusion is (1) Unionville granodiorite, a basic hypersthene-bearing granodiorite 
which itself has developed basic faces of granogabbro; (2) Clancy granodiorite; (3) 
porphyritic granodiorite; (4) biotite adamellite; and (5) muscovitic biotite granite, Alas- 
kite and aplite are abundant and were presumably (but not yet proved) developed most 
abundantly during the final stages of bathylithic consolidation, The order of emplacement 
of the successive intrusives is in the order.of increasing silicity. 

The Boulder bathylith and its satellitic stocks have exerted extensive contact meta- 
morphism, both purely thermal and pyrometasomatic. Most notably, the Helena dolomite 
has been transformed into aphanitic tremolite-diopside hornfels to a maximum distance 
of 10,000 feet from the edge of the bathylith. The highest rank of metamorphism attained 
is in sillimanite-cordierite-microperthite hornfels, remarkable rocks that have formed at 
widely separated localities, In places the magma has reacted with limestone xenoliths with 
the result that the xenolith is surrounded by an aureole of augite granodiorite in place of 
the normal hornblende-biotite granodiorite. In other places the evidence appears to de- 
mand that the magma in depth had dissolved limestone. By this syntexis alkalic rocks were 
generated that range from mildly alkalic, such as the Priests Pass leucomonzonite and 
the syenodiorite of the large stock northwest of Helena, to strongly alkalic, as represented 
by the nepheline shonkinite occurring east of Montana City. 


INTRODUCTION 


It is customary in a Presidential address to deal with some general prob- 
lem by way of critical review and synthesis, but I will beg your indulgence 
this evening if I break with precedent. I propose to discuss some of the prob- 
lems that I have met in field work on the Boulder bathylith of Montana during 
the past ten years. In this work I have been indefatigably assisted by Mrs. 
Knopf, partly in reciprocity I think for having accompanied her in her field 
work in the Taconic region of eastern New York, which has been a geologic 
battle ground for nearly 130 years. The Montana work has been done on my 
own time and expense, except that some of the chemical analyses have been 
paid for from the research fund of my Sterling Professorship in Yale Univer- 
* Address of the retiring President of the Peninsula Geological Society, Palo Alto, Cali- 
fornia, May 3, 1956. 
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sity. Consequently, whatever I am to say this evening has needed only to have 
the approval of myself. 

The Boulder bathylith of Montana is a large plutonic mass of 1100 square 
miles area that extends 70 miles along the Continental Divide from Mullan 
Pass on the north to a point 20 miles south of Buite. It was named in 1897 by 
W. H. Weed of the U. S. Geological Survey, because the region it occupies is 
drained by the Boulder River and tributaries and because Boulder, the county 
seat of Jefferson County, is near the center of the bathylith. 

The first geologist to have some broad ideas about the great intrusive 
mass that later was named the Boulder bathylith was Waldemar Lindgren. In 
1883 Lindgren was a member of the Northern Transcontinental Survey, which 
was sponsored by the Northern Pacific Railway. In the course of his work 
Lindgren examined the northern end of the bathylith near Mullan Pass, which 
is where the Northern Pacific crosses the Continental Divide. He (Lindgren, 
1886) briefly described the granitic rock and recognized that it resembled the 
Jurassic granite of Nevada as described by Zirkel, the petrographer of the 
Fortieth Parallel Survey. Furthermore, Lindgren found that coal-bearing beds 
considered to be of Laramie age rest on the eroded surface of the granite at 
Mullan Pass. He concuded therefore that the granite was possibly of Jurassic 
age. Later, in 1891, he conceded that “The possibility is not excluded however 
that the granite is of very late Mesozoic age”. 

Possibly at this time there began to develop in Lindgren’s mind a bold 
idea indeed, an idea that has strongly influenced the thinking of many geol- 
ogists. In 1915 Lindgren expressed the idea thus: “the fundamental fact in 
the Cordilleran region is that the igneous activity began along the present 
Pacific coast line and gradually extended eastward”. Twenty years later in 
referring to the Eocene intrusions of the Black Hills in South Dakota, he 
(1933, p. 158) said that “it is thought that they hail from the same magma 
which furnished the Coast batholith, its activities gradually spreading east- 
ward”, 

The first modern geologic work on the bathylith was done by the U. S. 
Geological Survey near the close of the nineteenth century in connection with 
studies of the mining districts by S. F. Emmons and especially by W. H. Weed 
and his assistants, notably Joseph Barrell, at Butte, Elkhorn, Helena, and 
Marysville. It was soon discovered that the youngest rocks invaded by the 
bathylith consist of a series of andesitic lava flows and breccias, which form 
its cover or roof. The andesites were thought to be Eocene and therefore the 
bathylith was considered to be Miocene. In the first flush of enthusiasm the 
exciting fact was announced that the intrusive mass had reached within a 
thousand feet of the Earth’s surface (Weed, 1902, p. 452). But greater famil- 
iarity with the region destroyed this early confidence and less definite estimates 
were later reached. 

Weed (1899) termed the rock that makes up the bathylith a quartz 
monzonite, which he regarded as being of very uniform composition through- 
out its entire area, except at contacts especially with limestone where it is 
basic. He regarded the bathylith as the result of an intrusion of a single 
magma, which has locally developed facies. In his final report on the Butte 
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district, Weed (1912) developed these ideas in somewhat greater detail. These 
ideas remained the general concept of the Boulder bathylith until the present 
investigation, which has shown that, like most other great bathyliths, it is a 
composite mass made up of successive intrusions of differing composition. 

In 1911 Dr. Lindgren, who was then Chief of the Metalliferous section 
and became Chief Geologist of the U. S. Geological Survey in that year, as- 
signed me to a field project in Montana; namely, to make a map of the 
Boulder bathylith and surrounding area from Marysville on the north to Butte 
on the south and to prepare a report on the ten mining districts in the region. 
I was given also a field assistant, Henry G. Ferguson, something I had never 
had during my five years service on the Alaskan Division; and to have at- 
tained to this high privilege gave a great lift to the morale of a young aspir- 
ing geologist. At any rate the report was completed early in 1912 and was 
issued by the Survey in 1913 (Knopf, 1913b). 

In retrospect the field season of 1911 has always seemed particularly in- 
teresting for a number of reasons, not the least of which was the exhilaration 
of field work in the summertime of Montana. In 1938 I decided to examine 
the Boulder bathylith as a possible problem in the mechanics of intrusion. A 
rapid reconnaissance indicated that its granitic rocks show no density strati- 
fication; for example, granodiorite from the top of Colorado Mt. at 7400 feet 
altitude has a density of 2.71 and similar granodiorite from a depth of 3600 
feet in the Stewart mine at Butte has a density of 2.72. It became apparent 
that in order to begin to understand the bathylith it would be necessary to 
map it petrographically, which had never been done, and furthermore to de- 
termine the stratigraphy of the surrounding strata, their structure, and their 
contact-metamorphic alteration. The ingrained belief that the Boulder bathy- 
lith is a quartz monzonite mass and that therefore all the rock that makes it 
up and all its satellitic stocks is quartz monzonite has led to the result that on 
extant maps gabbro, granogabbro, tonalite, granodiorite, adamellite, granite, 
graphophyre, alaskite, leucomonzonite, syenodiorite, and hornfelsed diorite 
porphyry are all shown as “quartz monzonite”. To distinguish these different 
rocks and to map them is a difficult and time-consuming task and has often 
proved to be highly frustrating, but if our understanding of the Boulder 
bathylith is to be increased the task must be done. 

In carrying out this mapping a specific gravity balance by which the 
density of hand specimens could rapidly be determined proved to be indis- 
pensable. The specimens collected each day were weighed and their densities 
determined, and the results served as valuable checks on the field determina- 
tions, This was true not only for the igneous rocks but also for the contact- 
metamorphic and the sedimentary rocks. 


THE STRATIGRAPHIC COLUMN 
The rocks invaded by the Boulder bathylith range in age from Precam- 
brian to late Cretaceous. The oldest beds are of Beltian age, the Helena region 


containing many of the type localities of the Belt series as named by Walcott 
(1899) in his pioneer studies. 
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The oldest Beltian unit in the Helena region is the Spokane formation 
comprising argillite and argillaceous siltstone, with sandstone interbeds, gen- 
erally of some pronounced red color—maroon, brick-red, Indian red, or lav- 
ender rose. Presumably, but not surely, four Beltian units occur below the 
Spokane formation, as the Helena region is 60 miles west of where the Neihart 
quartzite, the lowermost Beltian formation, crops out. 

Above the Spokane is the Empire formation, also named by Walcott; it 
is 1000 feet thick and consists of pale-green argillites alternating with mud- 
cracked rose or lavender argillites. 

Overlying the Empire formation is the Helena dolomite (Knopf, 1950, 
p. 837-838), originally named Helena limestone by Walcott. It consists mainly 
of 4000 feet of buff-weathering siliceous dolomite. A chemical analysis of a 
representative sample of the dolomite is given in table 1. Limestone and 
mixed carbonate rocks occur in minor amounts. Remarkable algal reefs up to 


TaBLe 
Analyses of Helena Dolomite and Other Rocks 


from the Helena Region, Montana 


I II Ill IV V 
SiO. 32.88 71.81 65.50 70.15 65.09 
Al.O; 3.48 14.63 16.84 14.41 14.20 
Fe.0; 16 82 1.62 1.68 87 
FeO 1.93 84 1.44 1.55 1.10 
MgO 12.89 37 a7 63 43 
CaO 18.42 1.21 1.60 2.15 5.00 
Na,O 73 4.02 4.58 3.65 2.60 
K.0 60 1.81 6.07 4.50 76 
H,O 72 33 20 68 6.38 
H.O 09 04 03 2.36 
TiO, 14 37 ot 42 47 
P.O; 05 05 18 12 08 
CO, 27.7 24 07 05 
MnO 10 08 ll 06 04 
SrO 02 06 
BaO 13 AS 17 
S 03 
.05 
99,90 99.77 99.82 100.00 99.65 
LessO=S LessO—F .02 
99.81 99.63 
d 2.78 2.63 2.65 2.38 


I Helena dolomite; west of B. M. 4450, Scratchgravel Hills, Analyst, E. H. Oslund. 
II Dellenite; north of Moose Creek, Tenmile Creek. Analyst, J. J. Engel. 


III Dellenite vitrophyre welded tuff; Continental Divide, north of MacDonald Pass. 
Analyst, Doris Thaemlitz. 


[V Average dellenite; compiled by Nockolds (1954). 


V_ Obsidian-like welded tuff; Great Northern Railway quarry, Wolf Creek. Analyst, 
E. H. Oslund. 
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15 feet thick and traceable for thousands of feet along the strike occur from 
the bottom to the top of the Helena dolomite. See plate 1. Algal growths of 
e this kind were originally called cryptozoons by Walcott, but in 1914 he re- 
named the genus Collenia, after M. Collen, of White Sulphur Springs, Mon- 


PLATE 1 
Collenia, in the Helena dolomite. 


A. Third biostrome from top in a succession of four biostromes occurring through 
a stratigraphic thickness of 75 feet. North side of Cherry Valley. northwest of Helena. 


¥ B. Largest Collenia so far found in the Helena dolomite. Bedding dips 15°E, West 
side of ridge on west side of Broadwater Gulch, 2 miles west of Helena. 
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tana, who had supplied him with much material. They were reported to be 
abundant in the Newland limestone (considered to be stratigraphically far 
below the Helena dolomite) from base to top, and to occur in the Spokane 
“shales”, but were not mentioned as occurring in the Helena dolomite. The 
presence of Collenia reefs or biostromes throughout the entire 4000 feet of the 
Helena dolomite indicates that the beds accumulated in a shallow basin never 
deeper than the depth at which photos, nthetic activity was possible, 

Overlying the Helena dolomite is the Marsh formation, 3000 feet thick 
(Knopf, 1950, p. 838). It consists of deep-red argillites, mud-flake breccias, 
and quartzites. Casts of salt crystals occur locally, and this is probably the 
one distinctive feature of the Marsh formation. Above the Marsh, forming 
the topmost member of the Belt series in this region is a quartzite belt 1800 
feet thick, which has been named the Greenhorn Mountain quartzite (Knopf, 
1950, p. 839), from the great mountain forming the Continental Divide north 
of Mullan Pass. 

The Belt series in this region aggregates 12,000 feet in thickness and is 
unmetamorphosed except near the Boulder bathylith. 

The Paleozoic begins with the Flathead quartzite of Middle Cambrian 
age. Between the Flathead and the Precambrian there is widespread discon- 
formity. At the Continental Divide the Flathead quartzite rests on the Green- 
horn Mountain quartzite; eastward at Helena it rests on Helena dolomite; still 
farther east, in the Belt Mountains, the Helena dolomite has been thinned to 
a feather edge. At least 5000 feet of strata were removed by erosion in the 
Helena region before the Flathead quartzite was laid down. Deiss (1940, 
p. 1092) thinks that 20,000 feet were removed, his argument being in part 
that mild orogenic activity and epeirogenic uplift took place in post-Beltian 
time, whereby mountains, which he called the Helena Mountains, were formed. 
These mountains according to Deiss were reduced to a peneplane before the 
deposition of the Flathead quartzite in Middle Cambrian time. The important 
evidence adduced by Deiss that extended erosion took place between the end 
of Beltian deposition and the beginning of Paleozoic sedimentation is one of 
several powerful reasons for not accepting the statement often made in regard 
to the Beltian geosyncline that a geosyncline can subside to a depth of 70,000 
feet before it is subjected to orogenic appression. 

The post-Beltian section, extending from Middle Cambrian to late Cre- 
taceous, is 9000 feet thick. It is not the result of continuous sedimentation, 
however, for at least two long interruptions occurred. The first break lasted 
from the end of the Cambrian to late Devonian time. The second interruption 
lasted from late Permian to middle Jurassic time. During the second interval 
the region was reduced to a peneplane, according to Condit (1918) and Collier 
(1927). 

Marine sedimentation began again late in Middle Jurassic time, signified 
by the deposition of the Ellis group. This event is held to mark the beginning 
of a new geosyncline, originally called Coloradoan by Schuchert (1939) and 
later renamed Rocky Mountain geosyncline by Dunbar (1949). 

In Colorado time rhyolite tuff was erupted, which marked the beginning 
of volcanism in this part of the Rocky Mountain geosyncline (Lyons, p. 451). 
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Later in the Cretaceous a thick series of andesites and latites was poured 
out. Eruption took place under submarine conditions, as shown by marine 
fossils found in the associated tuffs on Indian Creek on the east flank of the 
Elkhorn Mountains. These fossils were kindly determined for me in 1947 by 
Dr. J. B. Reeside, who identified them as of Niobrara age. 

Volcanism was later renewed subaerially. During this second episode 
andesites were erupted as well as abundant highly flow-banded rocks that 
resemble sparsely porphyritic biotite andesites, “pheno-andesites”. They carry 
scattered phenocrysts of andesine (Any) and a few of biotite. They cor- 
respond to dellenite, as shown in table 1. Many of them are welded tuffs. A 
chemical analysis of a black vitrophyric welded tuff of dull pitchy luster from 
the Continental Divide north of MacDonald Pass is given in table 1 (analysis 
III). It proves to be unusually high in K,O0—6.07 per cent. It corresponds 
closely in chemical composition to the silicic extrusive named toscanite by 
Washington (1897, p. 37), but called rhyodacite, quartz latite, dellenite, 
subaluminous rhyolite, as well as other names by other petrographers. Its 
chemical composition corresponds to that of an orthoclase-rich quartz latite 
but neither orthoclase nor quartz have developed modally. 

That welded tuffs were formed during the second period of volcanic acti- 
vity in late Cretaceous time proves that the eruptions of this period, in contrast 
to the earlier submarine outbreaks, took place under subaerial conditions. In 
the Adel Mountain volcanics, a series of potash-rich lavas occurring on the 
Missouri River north of the Boulder bathylith, | found brilliantly vitreous 
obsidian-like welded tuffs (Barksdale, 1951). An analysis is given in table 1 
(analysis V); it shows that the analyzed rock is of unusual composition and 
indicates that further investigation of the Adel volcanics is necessary. The 
obsidian-like welded tuffs are associated with green and purple lithoid welded 
tuffs. These volcanics had earlier been mapped by Pardee and Schrader 
(1933) as part of the “andesite-latite” series that consists of the youngest 
rocks invaded by the Boulder bathylith. This correlation is probably strength- 
ened by finding that welded tuffs occur both in the roof rocks of the bathylith 
and in the Adel Mountain volcanics. The correlation is of particular signifi- 
cance, as Lyons (1944, p. 460) found leaves in the Adel Mountain volcanics, 
from which their age is inferred to be very late Cretaceous, probably late 
Montanan. The major diastrophism, as emphasized by Lyons, took place after, 
presumably soon after, the eruption of the Adel Mountain volcanics. 


ONSET OF THE LARAMIDE OROGENY 


The main deformation that has affected the region now set in, presum- 
ably, as just mentioned, shortly after the final volcanism. It was the first severe 
Laramide orogenic crisis in Montana. The Beltian, Paleozoic, and Mesozoic 
rocks, including the youngest volcanics, were folded, involving a thickness of 
20,000 feet or more. These rocks had accumulated in three geosynclines, cur- 
rently named the Beltian, Cordilleran, and Rocky Mountain. The Beltian and 
Cordilleran geosynclines were of the kind called miogeosynclinal. The Rocky 
Mountain geosyncline, however, which began its life late in the Jurassic and 
ended late in the Cretaceous after the eruption of the two series of volcanic 
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rocks just mentioned, was manifestly pliomagmatic in this part of Montana 
and can probably qualify as a eugeosyncline. 

The regional folding was of an open character, with the resultant folds 
trending northwest. The folding was followed or accompanied by the intrusion 
of the Boulder bathylith, which was emplaced in successive pulses or surges 
of magma of differing compositions. 

Near the bathylith the invaded country rock has been more closely folded 
than at a distance from the contact. In places the strata adjacent to the bathy- 
lith stand vertically and have even been overturned. Locally a series of reverse 
faults has developed along the eastern border of the bathylith. The intrusive 
magma manifestly has made room for itself partly by crowding aside the 
enveloping rocks, by close appression of the beds, by overturning them, and 
by imbricate high-angle thrusting. 

To Stille (1940, p. 264) the intrusion of the Boulder bathylith was the 
result of “allochthonous synorogenic plutonism fed by the lateral migration of 
sialic magma.” The magma had migrated from the west in the direction of 
the eastward progression of folding. The feebler intensity of folding in the 
Rocky Mountains compared to the Nevadan folding farther west precluded 
the possibility of autochthonous synorogenic plutonism. So, says Stille, “I 
have come by way of magmatectonic considerations to a conclusion that 
Lindgren had formulated from a more mineralogic-petrographic viewpoint.” 


SATELLITIC STOCKS 


Some 20 stocks occur in proximity to the northern end of the Boulder 
bathylith. They comprise a wide range of petrographic types: olivine-ortho- 
clase gabbro, diorite, tonalite, adamellite, granodiorite, granite, alaskite, 
syenodiorite, and diorite porphyry. Every one of them presents this problem: 
(1) was it emplaced before the intrusion of the bathylith or (2) was it em- 
placed during the intrusion of the bathylith, and if so, with which surge of 
magma; or (3) was it emplaced long after the intrusion of the bathylith? 
Most of the stocks appear to be genetically related to the bathylith, but to have 
been emplaced at different stages in its history. The best known of the satel- 
litic stocks is the Marysville granodiorite stock. In an earlier paper I have 
discussed in detail the problem of determining that the Marysville stock is 
genetically allied to the Boulder bathylith (Knopf, 1950). The largest stock 
is the Scratchgravel Hills stock, 8 square miles in area, which forms the group 
of hills northwest of Helena, rising out of the great circular intermontane de- 
pression known as the Helena Valley. It proves to be an anorthoclase sye- 
nodiorite of weakly alkalic character. Its chemical composition is shown by 
analysis VII in table 2. 

Recently Rittman (1953) has shown that from the chemical analysis of a 
rock can be calculated an index “p”, which gives a measure of the alkalicity’ 
of the rock. The dividing point between lime-alkalic and alkalic rocks is at 55; 
values between 50 and 55 indicate a weakly alkalic character, and values be- 
tween 55 and 60 indicate a weakly lime-alkalic character. The value of the 


1 “Alkalicity” is suggested as a more appropriate term than “alkalinity”, which is a 
chemical term and is measured by pH. 
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TABLE 2 
Analyses of Plutonic Rocks of the Boulder Bathylith and Two Related Stocks 
I II Ill IV V VI VII Vill IX xX 

SiO, 61.14 61.40 54.63 65.49 66.14 71.28 6848 56.01 6413 51.95 
Al:Os 15.28 1541 1660 1449 15.69 1450 1493 17.84 16.77 12.58 
Fe.0; 1900 «63.32 2.11 1.93 1.04 1.32 3.56 196 63.12 
FeO 452 S76 “Si “200 1,19 2.07 3.05 169 3.91 
MgO 3.41 299 465 245 ~~ 1.60 .98 1,37 2.57 129 6.62 
CaO 5.51 486 7.93 429 3.71 2.45 2.61 7.39 2.72 13.00 
Na,O 2.78 2.69 248 280 3.50 3.16 3.01 4.01 3.92 2.85 
K:.0 3.34 412 2.110 3.66 3.54 4.29 4.59 3.33 5.77 4.61 
H,0 + 65 85 56 AZ 43 .60 61 39 
H,O- 09 18 14 05 05 13 12 24 19 .03 
TiO: 82 80 89 65 43 26 39 65 2 
P.O; .26 .26 42 09 42 23 46 
CO, .08 02 05 08 02 

MnO 12 .09 .16 .10 .06 07 19 05 14 
SrO 04 07 04 04 08 02 02 .10 02 .06 
BaO .08 12 05 13 ll Al 12 05 
01 02 04 03 O01 01 

99.65 99.94 99.74 

Les O = S 02 02 

Total 100.03 99.64 99.94 99.91 99.73 100.10 100.17 100.09 99.78 99.88 
d 2.78 2.76 284 2.71 2.70 2.65 2.67 2.460 «62666-2393 


I Unionville granodiorite; type locality, old quarry 0.9 mi. NW of BM4589, Clark 
Gulch. Analyst, J. J. Engel. 


If Unionville granodiorite; near Benson xenolith. Analyst, Doris Thaemlitz. 
III Granogabbro; at contact NW of Park City. Analyst, Doris Thaemlitz. 
IV Clancy granodiorite; Kain quarry, Clancy Creek, Analyst, J. J. Engel. 


V_ Porphyritic granodiorite; Schimpf quarry, 7 mi. W of Helena. Analyst, Doris 
Thaemlitz. 


VI Biotite adamellite; ridge between:Jackson Creek and Lump Gulch. Analyst, J. J. 
Engel. 


VII Biotite adamellite; Broadwater stock, Helena, Analysts, B. Smith and R. B. Ellestad. 


VIII Syenodiorite; old quarry near Great Northern Railway, Scratchgravel Hills, Ana- 
lyst, L. C, Peck. SrO determined by S. S. Goldich. 


IX Priests Pass leucomonzonite ; road cut just west of Priests Pass. Anal st H. 
yst, 
Baadsgaard. 


X Nepheline shonkinite; 3/4 mi. southeast of Montana City. Analyst, H. Baadsgaard. 


alkalicity index for the syenodiorite is 53.8. For the prevailing rocks of the 
Boulder bathylith the value averages around 62. 

This stock of syenodiorite is the first recorded occurrence of alkalic rocks 
in association with a large mass of ordinary lime-alkali rocks within the Rocky 
Mountains of Montana. East of the Rockies is the remarkable alkalic province 
of the Highwood Mountains made famous by Pirsson (1905) and by Larsen 
(1940). 

AGE OF INTRUSION OF THE BOULDER BATHYLITH 

The time of intrusion of the Boulder bathylith has been narrowed down 

somewhat in recent years. In the marine volcanic tuffs on Indian Creek fossils 
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were found, which as already mentioned, are of Niobrara age, i.e., of late 
Colorado age. By correlation of the welded tuffs in the roof of the bathylith 
with the welded tuffs in the Adel volcanics, a very late Montanan or possibly 
post-Montanan age is indicated for the youngest rocks invaded by the bathy- 
lith. 

In the southwestern Montana, according to Eardley (1951), the first 
main phase of the Laramide orogeny took place in latest Cretaceous or early 
Paleocene time. As a result of this orogeny the very coarse and thick Beaver- 
head conglomerate was laid down. The Beaverhead conglomerate makes up 
the greater portion of the Beaverhead formation, of exposed thickness of 
9700 feet, and according to Lowell and Klepper (1953) the middle portion 
of the formation is tentatively indicated on the basis of fossil fresh-water 
gastropods to be of early Paleocene age. 

By the beginning of the Oligocene the bathylith had been deeply un- 
covered by erosion, for at Pipestone Springs, east of Butte, beds carrying a 
well-defined vertebrate fauna of earliest Oligocene—Chadronian—age (Wood 
et al., 1941) rest on the surface of the granodiorite. So far as the geologic 
evidence goes, the bathylith was emplaced sometime between late Cretaceous 
and the beginning of the Oligocene. 

Five determinations of the absolute age of the Boulder bathylith have 
recently been reported by Chapman, Gottfried, and Waring (1955). They are 
based on alpha-counting and measurement of the lead content of zircon and 
monazite obtained from rocks of the Boulder bathylith. The reported ages 
range from 62 to 72 million years. Dr. J. H. Reynolds, of the Physics Depart- 
ment of the University of California at Berkeley and Dr. R. E. Folinsbee, of 
the University of Alberta, have kindly determined for me, by the potassium- 
argon method, the age of the granodiorite that makes up the larger part of the 
Boulder bathylith (Knopf, 1956). Dr. Reynolds reports “that the most prob- 
able age is 87 million years”. This potassium-argon age determination suggests 
that the Bouder bathylith was emplaced late in the Cretaceous; consequently 
during the first main phase of the Laramide orogeny. The time of intrusion 
cannot be placed closer in the absolute geologic time scale however, until 
agreement is reached on what strata mark the end of the Cretaceous and until 
we have a definitive determination of the absolute age of those strata, be they 
topmost Maestrichtian or topmost Danian. 


COMPOSITE NATURE OF THE BATHYLITH 


The Boulder bathylith has heretofore been considered to be a one-magma 
bathylith. As the present mapping proceeded it became clear, however, that 
the bathylith was built up by successive intrusions, which arrived in place 
according to the well-known order of increasing silica content. In short, the 
bathylith proves to be a composite mass that has been built up by five or more 
magmas arriving one after the other in what the evidence suggests was a 
geologically short time. 

Unionville granodiorite.—The earliest intrusion is a basic granodiorite. It 
is dark, heavy, and resembles a diorite. Neither quartz nor potassium feldspar 
are visible to the unaided eye, but as seen under the microscope enough po- 
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tassium feldspar and quartz (15 percent) are present to allow the rock to 
quality as a granodiorite. It is a hypersthene-bearing augite-biotite-hornblende 
granodiorite. A chemical analysis of the type rock is given in table 2. 

There are at least seven definitions of granodiorite: by Lindgren (1900), 
Iddings (1913), Niggli (1931), Johannsen (1932), Tréger (1935), Hatch, 
Wells, and Wells (1949), and Nockolds (1954), but none agrees fully with 
any other. Niggli and Johannsen agree in drawing the line between granite 
and granodiorite so that granites are distinguished by having potassium 
feldspar predominate over plagioclase and granodiorite by having sodic 
plagioclase predominate over potassium feldspar. This usage has the advantage 
of simplicity, and as pointed out by Johannsen (1932) makes the Brogger 
term quartz monzonite [= adamellite| superfluous. 

Nockolds (1954) in a notable paper giving newly computed average 
chemical compositions of igneous rocks presents a tabular classification of 
igneous rocks. He distinguishes adamellite from granodiorite as follows: 
adamellite has a potassium feldspar content between 60 percent and 40 percent 
of the total feldspar content, whereas granodiorite has a potassium feldspar 
content between 40 percent and 10 percent of the total feldspar. In accordance 
with Tréger and Shand, 10 percent is taken as the lower limit for the quartz 
content of granite, adamellite, granodiorite, and tonalite. 

These definitions appear to be logical; furthermore they embody the 
greatest measure of agreement so far reached and will therefore be used in 
describing the granodiorite and allied rocks of the Boulder bathylith. 

The basic granodiorite is well exposed south and southeast of Helena. 
Hence it is proposed to call it the Unionville granodiorite, from a small village 
south of Helena. Near the contact it develops a basic facies having the high 
density of 2.85 or more. Its chemical analysis is given in table 2 (no. III). 
Petrographically the facies proves to be a hypersthene granogabbro, grano- 
gabbro being the name given by Johannsen (1932, p. 367) to the orthoclase- 
bearing variety of quartz gabbro. At many places the granogabbro and 
granodiorite contain numerous large plates of poikilitic biotite 1 to 3 cm. in 
diameter, the result of biotitization. That is, the biotite continued to grow afier 
the other constituents had ceased growing, and has spread out at the expense 
of the plagioclase. This continued growth of the biotite is interpreted as a 
magmatic end-stage reaction or deuteric effect. It has produced very striking- 
looking rocks of coarsely poikilitic habit, much more “basic” in appearance 
than they really are. 

Clancy granodiorite.—The second intrusion is the Clancy granodiorite, 
so named from the Kain quarry on Clancy Creek. It is a coarse-grained rock 
in which quartz is conspicuous. Its chemical analysis is given in Table 2 (no. 
IV). The earlier, Unionville granodiorite had already solidified and was cut 
by dikes of the younger invading magma, as is well shown on the summit of 
Colorado Mountain (plate 2). The Clancy granodiorite is not perceptibly 
chilled against the earlier granodiorite. However, the light-gray color of the 
younger, Clancy granodiorite contrasts strongly with the dark gray of the 
earlier, more basic, granodiorite. The Clancy granodiorite is in appearance, 
density, and chemical composition like the bathylithic rock at Butte, the so- 
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Clancy granodiorite (light gray) cutting the Unionville granodiorite, which is the 
earliest intrusive of the Boulder bathylith. Summit of Colorado Mountain. 


called Butte quartz monzonite. The naming of the Butte rock has had a re- 
markable history. Originally named the Butte granite by S. F. Emmons 
(1888), it continued to be called so by Weed in the Butte Special Folio (1897), 
for the reason that “thin sections of the rock studied under the microscope 
show that orthoclase is the most abundant mineral, and that the rock is there- 
fore a granite.” At the same time however, Weed recognized that the chemical 
analyses, of which four were available, “closely resemble those of the granod- 
iorites of the Gold Belt of California”. Later, Weed (1899) renamed the 
“granite” of Butte “quartz monzonite,” without much explanation. 
Subsequently detailed work was done by Cross, Iddings, Pirsson, and 
Washington (1903) on the microscopical petrography of the Butte “granite”. 
Two thin sections of the granite, one from rock near the Parrott vein and the 
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other of unknown provenance, were measured according to the Rosiwal 
method, with closely agreeing results. They found that the zoned plagioclase, 
considered to average Ano in composition, amounted to 41.98 percent by 
weight and the orthoclase to 18.7 percent. Consequently, if they had followed 
the proposal of Lindgren (1900) to redefine granodiorite as having sodic 
plagioclase characteristically “at least double the amount of the alkali feld- 
spar,” they would have called it granodiorite ; however, they remained content 
to refer to it as “granite”. They computed also the mineral composition of the 
“granite” at Walkerville, on the north edge of Butte, which as well as its horn- 
blende and biotite had been chemically analyzed. This computation also 
showed that the amount of plagioclase (Any) is more than double the ortho- 
clase. 

Despite this demonstration, Weed (1912, p. 32) in his final memoir on 
the Butte district reaffirmed that the rock at Butte is a quartz monzonite, be- 
cause “thin sections of the rock studied under the microscope show that the 
orthoclase and plagioclase feldspars are present in nearly equal amounts, so 
that the rock according to precise petrographic nomenclature shou!d be called 
a quartz monzonite”. 

An attack from another angle on the use of the term quartz monzonite 
comes from a detailed restudy of the rocks of the type area Predazzo-Monzoni, 
in the Tyrol (Leitmeier, 1940). He affirms that Brégger ought not to have 
proposed to use “quartz monzonite” for rocks other than those called quartz 
monzonite by Romberg and by Rosenbusch, which are rocks that contain not 
over 5 percent of quartz. All the quartz monzonites in Predazzo-Monzoni ac- 
cording to Leitmeier (1940, p. 180; 1948, 79-81) “are true monzonites and 
have no kinship with the quartz-rich alkalicalcic granitic rocks of Walkerville, 
Butte, Montana, erroneously called quartz monzonite”. 

Porphyritic granodiorite—The third intrusion is here called the por- 
phyritic granodiorite. It is characterized by its content of orthoclase pheno- 
crysts an inch in length and 1/2 inch across. It is more constant in composition 
(table 2, no. V) than any of the other intrusive masses that make up the 
bathylith. Its average density is 2.68, the range being only from 2.71 to 2.65. 
The plagioclase is zoned from Ang; in the core to Ano in the peripheral zone. 
Granodiorite of this kind extends to Mullan Pass, the northernmost limit of 
the bathylith. 

The position of the porphyritic granodiorite as the third member in the 
intrusive sequence remains so far inferential, as critical contacts have not been 
found. 

Biotite adamellite——The fourth intrusion is a “white granite”, as it was 
called in the field; petrographically, however, it is a biotite adamellite. It is 
definitely intrusive into the Clancy granodiorite, and a dike 5 feet thick was 
found to cut the Clancy granodiorite. The dike has unchilled borders, so 
manifestly no great interval of time separated the successive intrusions. 

In chemical composition the adamellite (table 2, no. VI) is almost 
identical with that of the average biotite adamellite as recently computed by 
Nockolds (1954). The adamellite is higher in silica (71.3 percent), much 
higher in quartz, and is lower in density (2.65) than the earlier intrusives. 
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The plagioclase is Any however, like that in the porphyritic granodiorite. 
Therefore the evolution of the bathylith as measured by the anorthite content 
of the plagioclase did not advance during the interval between the third and 
fourth intrusives. In terms of the ferromagnesian minerals, however, the evolu- 
tion had advanced so far that only biotite was crystallizing from the magma. 
The type locality is on the ridge between Jackson Creek and Lump Gulch. 
Large masses occur on Clancy Creek; on U. S. Highway 91 south of the 
Helena-Boulder divide; and along U. S. 10S from one mile east of Pipestone 
Pass to the water fountain on the west side of the Continental Divide. 

Biotite granite—The fifth intrusive is a medium-grained biotite granite 
low in plagioclase (Anzo). Muscovite is generally present, but it is all of re- 
placement origin and is probably deuteric. This granite (d = 2.61) forms 
relatively small masses and clearly intrudes the Clancy granodiorite, into 
which it sends dikes. It is in places remarkably and conspicuously tourmalini- 
ferous, being crowded with tourmaline nodules, each nodule being surrounded 
by a white border 1 cm. wide. These features are spectacularly developed near 
Montana City southeast of Helena. Younger than the muscovite-biotite granite 
are aplite dikes up to 40 feet thick, cutting both the granite and the country 
rock. The aplites are highly tourmaliniferous, the tourmaline occurring both 
in disseminated form and localized along fractures. 

The muscovite-biotite granite was mapped by previous investigators as 
aplite. The discovery that it is itself cut by dikes of genuine aplite (d = 2.59) 
is therefore significant and proves that these aplites represent a still more 
advanced stage in the evolution of the bathylith. 

Alaskite and aplite (aplite might well be called microalaskite) are com- 
mon throughout the area of the Boulder bathylith. In fact, aplite has been 
estimated by one investigator (Knopf, 1913b, p. 34) to make up 5 percent of 
the exposed area of the bathylith and 10 percent by another investigator 
(Billingsley, 1915) ; but these estimates were not supported by areal mapping 
or by precise petrographic knowledge. Alaskite and aplite have doubtless been 
formed during each of the five recognized stages of the development of the 
bathylith, but they appear to have been more abundantly developed in con- 
nection with the later, more silicic intrusives. Inasmuch as an aplite, regard- 
less of which of the five intrusive masses to which it is genetically linked, 
represents an ultimate stage in differentiation, it was of interest to compute 
the new alkalicity index for the granodiorite at Walkerville, Butte, containing 
63.88 percent of SiO., and for the aplite (76.87 percent of silica) cutting it. 
For both these rocks excellent chemical analyses are available. The indices 
prove to be almost identical (61.1 for the granodiorite and 60.4 for the aplite), 
and indicate therefore that the extreme differentiation necessary to produce 
the aplite has not changed the value of the alkalicity index. 

ABERRANT INTRUSIONS 

The five intrusions, ranging from Unionville granodiorite to muscovite 
granite, constitute the normal stages of the evolving bathylith. Some other 
masses represent unsolved problems. 

One such mass occurs.at Priests Pass just south of Mullan Pass, and is 
unlike the rocks already described. It is of monzonitic composition nearly 
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devoid of quartz and has but a very moderate content of biotite and horn- 
blende. It is here called “Priests Pass leucomonzonite.” Its chemical composi- 
tion is given in table 2 (no. IX). The contact between the leucomonzonite 
and the porphyritic granodiorite is cleanly and perfectly exposed southwest 
of Austin, and gives convincing evidence that the leucomonzonite is not a 
facies of the granodiorite. Inasmuch as neither rock is chilled or otherwise 
modified against the other, the time sequence remains unknown. 

Another mass but different from the Priests Pass leucomonzonite, con- 
stitutes the high mountainous portion of the bathylith south and east of 
Helena. It consists of a darkish granitic rock resembling to the eye the Union- 
ville granodiorite. Unexpectedly it proved to be, as seen under the microscope, 
a graphophyre, or “granophyre”, with phenocrysts of labradorite and py- 
roxene in a superbly developed groundmass of micropegmatite. In reference 
to “granophyre”, Michel-Levy (1889, p. 19) said that the name “me parait 
aussi mal choisi que possible pour représenter des roches dont le magma 
présente les formes de la pegmatite graphique.” No other rock mass of this 
kind has previously been reported from the Boulder bathylith. 

The position of the graphophyre in the intrusive sequence is only partly 
known: it cuts the earliest intrusive, the Unionville granodiorite, in dikes and 
is itself cut by tourmaliniferous aplite. Mapping the graphophyre has proved 
to be an exacting and nearly impossible task. 


CONTACT METAMORPHISM 


The metamorphism exerted by the Boulder bathylith on its enveloping 
rocks is a notable feature. According to Billingsley however, (1915, p. 47), 
“the intrusion produced only slight metamorphic effects”, and Grout and Balk 
(1934) wrote “For a coarse-textured intrusive with an area of about 1200 
square miles, the zone of prominent contact metamorphism is remarkably 
narrow—at some places it is hardly recognizable”. To which the answer can 
be made “Never underestimate the power of a bathylith”. 

The voluminous development of calc-hornfelses in the contact zones is 
impressive. The calc-hornfelses occur in nine formations at least, ranging from 
the Precambrian Empire formation to the Cretaceous. They are white, light- 
gray, or creamy-colored aphanitic rocks of chert-like aspect. They consist of 
diopside and tremolite, and subordinate scapolite. The Helena dolomite has 
furnished the largest amount of calc-hornfels, extending as far as 10,000 feet 
from the visible contact. Despite the drastic metamorphism, the weathered 
outcrops of hornfelsed Helena dolomite retain completely the innocent ap- 
pearance of normal sedimentary rocks (plate 3). 

An intensive study was made of the contact zone that surrounds the 
Broadwater stock in the western outskirts of Helena. Here the Helena dolomite 
has been transformed into an aureole of aphanitic calc-hornfels (or better 
calc-magnesian hornfelses) that encircles the stock and ranges from 300 to 
2500 feet in width. The chief minerals are diopside and tremolite, with wollas- 
tonite occurring only at the immediate igneous contact. In the outer fringe 
of the aureole sodic scapolite and the white chlorite, leuchtenbergite, have de- 
veloped. A bed of wholly unmetamorphosed blue odlitic limestone occurs well 
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PLATE 3 


Cale hornfelses (hornfelsed Helena dolomite) in the contact-metamorphic aureole 
of the Broadwater stock, showing retention of bedding and sedimentary appearance. Four 
miles west of Helena, on U. S. Highway 10N. Photograph by Geo. Switzer. 


within the contact zone, where it is intercalated between rocks completely 
altered to calc-hornfels; this intimate association of unaltered and completely 
altered rocks emphasizes vividly that dolomite is enormously more susceptible 
to thermal metamorphism than limestone. 

The widespread development of these calc-hornfelses by the Boulder 
bathylith and its satellitic stocks strikingly illustrates the principle of dedolo- 
mitization announced by Teall (1903)—namely that during metamorphism of 
dolomite the magnesium is liberated as periclase (MgO), or if silica is avail- 
able the magnesium is preferentially silicated, resulting in the liberation of 
calcium carbonate. 
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Cordierite hornfels is another notable feature of the contact metamorph- 
ism. No less than seven formations have developed cordierite hornfels. The 
most striking development is shown by the Park argillite, 200 feet thick, which 
mile after mile has been converted solidly into cordierite hornfels. A show- 
place is east of Austin on the Mullan Pass highway. Large blocky outcrops of 
massive black rock resembling a basic igneous rock occur here. On lightly 
weathered surfaces of the hornfels the cordierite resembles wheat grains or 
fusulines, but on strongly weathered surfaces the cordierites have been etched 
out and the hornfels closely simulates a vesicular basalt. Among the many 
geologists to whom | have shown this cordierite hornfels was Prof. Goodspeed, 
of the University of Washington, who took back to the University a cubic foot 
or so. On the field trip during the Butte meeting of the Rocky Mountain 
section of the Geological Society of America in 1953 this show place was one 
of the stops. The geologists present were much impressed by the cordierite 
hornfels, but a group of graduate students from the University of Washington 
were particularly interested: they saw at last what the rock they had been 
calling vesciular basalt in their PhD orals really was and where it came from! 

The highest rank of metamorphism attained by any of the rocks in the 
contact zone is represented by sillimanite-microperthite-cordierite hornfelses. 
These are notable rocks in which the sillimanite occurs in minute gleaming 
}risms in random arrangement, and the microperthite is distinguished for its 
regular, evenly distributed, delicate exsolved discs of albite. Probably on 
further study this microperthite will serve as a better geologic thermometer 
than the sillimanite. These rocks occur only in roof pendants or in immediate 
contact with the bathylith, and although they have been found at widely 
separated localities they occur only near the contacts with the Unionville 
granodiorite, the earliest and presumably the hottest intrusion. 


REACTION OF THE BATHYLITH WITH THE COUNTRY ROCKS 


On my preliminary examination of the bathylith I found a large perfect 
exposure of a portion of its extreme northern contact. A broad apophysis of 
the porphyritic granodiorite extends conveniently across the Mullan Pass Road 
into the Madison limestone. The granodiorite in this apophysis, even at the 
borders and tip of the apophysis, is completely normal, its density is un- 
changed, its orthoclase phenocrysts are present in normal fashion—manifestly, 
the magma absorbed no limestone, nor has it been basified or otherwise 
changed. The Madison limestone has been transformed to a diopside-bearing 
marble, which also is the normal behavior: diopside at the immediate contact 
with the bathylith, tremolite at a greater distance. Only later did it become 
plain that the magma has in other places reacted with limestone, with results 
of a remarkable character and of great diversity. “Never underestimate the 
power of a bathylith”—this was learned the hard way. 

At five widely separated localities magma and limestone have reacted 
with striking results. 

The first example is at the Benson ranch, southeast of Helena. At this 
locality a xenolith, 40 feet long and 20 feet wide at most, is enclosed in the 
Unionville hornblende-biotite granodiorite. As shown by its chemical analysis 
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in table 2 (no. II), the granodiorite has closely the composition of the Union- 
ville granodiorite at the type locality; its density is 2.76 and its plagioclase 
has the composition An,,. The xenolith is made up mainly of andradite, py- 
roxene, epidote, and minor highly calcic plagioclase (Ansgo-90). The granodio- 
rite that surrounds the xenolith is lighter-colored than the normal granodiorite 
to a distance of 40 feet from the xenolith and is characterized by numerous 
patches of poikilitic orthoclase. Unexpectedly, this lighter-colored granodiorite 
has the high density of 2.90 next to the xenolith and shades off farther away 
to 2.80. The increased density is due to a high content of augite in place of 
hornblende and biotite. These relations are interpreted thus: the xenolith 
originally consisted of limestone or dolomite; part of it was dissolved by the 
magma; and the resulting increase of CaO and MgO in the magma suppressed 
the formation of hornblende and biotite and caused augite to form in their 
stead, and the potash that would have gone to make biotite was used to form 
the poikilitic orthoclase. In short, an augite granodiorite was formed instead 
of a hornblende-biotite granodiorite. Remarkably, the augite granodiorite in 
immediate contact with the tactite has a plagioclase of composition only as 
calcic as Anso, whereas the plagioclase in the tactite has the composition Ango. 
Finally, after the augite granodiorite had consolidated, pyrometasomatic 
activity ensued and converted the residual portion of the xenolith into a tactite 
composed of clinopyroxene (aegirine-augite), andradite, bytownite, and 
epidote. Minor actinolite, quartz, calcite, magnetite, and analcite and traces of 
pyrite also occur. The andradite has replaced clinopyroxene, and epidote has 
replaced both clinopyroxene and andradite, and has spread in very minor 
amount into the adjacent augite granodiorite. 

Apparently the effect of the assimilation of limestone by the granodiorite 
magma at the Benson xenolith is somewhat analogous to phenomena at Tenn- 
berg, Sweden, described by von Eckermann (1923). There the granite ad- 
jacent to a limestone xenolith was rendered non-biotitic for a distance of 200 
feet from the xenolith. But instead of augite having formed as the result of 
the addition of lime to the magma as in Montana, the hastingsite variety of 
hornblende was formed. Tilley (1949) also has found that at dolomite skarns 
the hornblende and biotite of the adjacent granite have given place to clino- 
pyroxene. The suppression of biotite is probably the result of the strong af- 
finity of lime for alumina to form the anorthite molecule; as long as any free 
lime was present in the magma biotite could not form. 

More spectacular effects than at the Benson xenolith were produced at 
a xenolith or roof pendant east of Montana City. This xenolith consists of 
limestone and is 150 feet long by 30 feet wide at its widest portion. Attention 
was <irected to this xenolith by the wollastonite and coarse yellow vesuvianite 
that form large masses gorgeously colored in the brilliant light of a Montana 
summer. Less conspicuous associates are grossularite and diopside. One end 
of the xenolith, the south end, has been converted solidly into pyrometaso- 
matic rock, whereas the north end still consists of bedded limestone. The 
pyrometasomatic rock is traversed by dikes of nepheline shonkinite, plagio- 
clase shonkinite, and allied rocks. 
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Chemical analysis of the nepheline shonkinite (table 2, no. X) shows 
that it is like that of the type shonkinite from Square Butte in the Highwood 
Mountains and that at Shonkin Sag, Montana. The alkalicity index of the 
nepheline shonkinite is 42, confirming, if that were necessary, its highly al- 
kalic character. Here then is powerful support for Daly’s theory that alkalic 
rocks are generated from lime-alkalic magmas by the absorption of limestone. 
You will recall that Dr. Tilley (1931) found at Scawt Hill, Ireland, that 
nepheline-bearing rock had been formed in a contact zone by reaction of a 
dolerite dike cutting chalk, the first field demonstration of the probability of 
Daly’s theory of the origin of alkalic rocks by absorption of limestone. When 
[ showed Dr. Tilley the Montana occurrence he exclaimed “you have found 
the underground equivalent of Scawt Hill”. It must be admitted, however, that 
the place where the nepheline shonkinite magma was formed is still deeper 
underground. 

This formation of shonkinite east of Montana City as a reaction product 
between granodiorite magma and carbonate rock is not a unique example. On 
the border of a small granodiorite stock on Grizzly Gulch just southwest of 
Helena there are shonkinite dikes up to 2 feet thick that cut a mass of vesu- 
vianite-diopside tactite 100 feet thick and several hundred feet long. 

The most voluminous production of shonkinite took place near Austin, 
at the northern border of the Boulder bathylith. The Madison limestone here 
has been converted into diopsidic tactite through a thickness of several hun- 
dred feet. Intrusive into this tactite are bodies of shonkinite as much as three 
hundred feet thick. The shonkinite in turn is cut by thin dikes composed 
wholly of coarse-grained alkali-feldspar rock—leucosyenite, or perthosite 
(Phemister, 1926). 

In all the places where shonkinite has so far been found near the Boulder 
bathylith this triple association occurs: limestone, tactite, and granodiorite. 


TABULAR SUMMARY OF BOULDER BATHYLITH PHANERITES 
AND ASSOCIATED ALKALIC ROCKS 

Table 3 gives a summary view of the evolution of the plutonic rocks of 
the Boulder bathylith. It begins with the Unionville granodiorite and its basic 
border facies the granogabbro from the contact at Park City southwest of 
Helena, and ends with the muscovitic biotite granite at Montana City. The 
table brings out clearly the strong contrast between the normal rocks of the 
bathylith and the alkalic rocks: the syenodiorite of the Scratchgravel Hills 
stock, the Priests Pass leucomonzonite, and the nepheline shonkinite, which 
are regarded as having been formed by limestone syntexis. 


LOCALIZATION OF ORES IN THE BOULDER BATHYLITH 
A major problem that the Boulder bathylith presents is what are the 
factors that have determined the remarkable localization of the ore deposits. 
The supreme outstanding feature of the ore deposits that are genetically re- 
lated to the Boulder bathylith is the enormous concentration of copper and 
zinc at Butte. Another facet of this problem is what factor or factors have 
determined the localized or restricted distribution of certain metals to certain 
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TABLE 3 
Parameters of the Boulder Bathylith Phanerites and Associated Alkalic Rocks 


SiO Alkalicity An Normative Color 
content Density index in plag. plag. qtz. index 
l nionville 
Granodiorite 61.14 2.78 60.7 17 44 15 21.6 
1A Granogabbro 54.63 2 84 62.2 55 56 8.6 28.9 
2 Clancy . 
Granodiorite 65.49 2.71 62.1 44 39 22 15.0 
Porphyritic 
Granodiorite 66.14 2.70 61.6 10-20 34 22 10.4 
t Adamellite 71.28 2.65 62.8 410 29 33 5.8 
5 Granite 2.61 30 
Syenodiorite 56.01 2.76 53.8 35 31 3.0 21.8 
Priests Pass 
Leucomonzonite 64.13 2.66 53.2 37 24 11.8 8.7 
Nepheline 
Shonkinite 51.95 2.93 12.3 48.1 


areas in the bathylith. Outside of Butte the ore deposits have been mainly 
silver-bearing lead, in many places associated with abundant tourmaline, com- 
prising the tourmalinic silver-lead type of ore deposit (Knopf, 1913a). Yet 
in other places in the bathylith extraordinary deposits of tourmaline occur, 
but no ore deposits were formed. On the other hand tourmaline is practically 
unknown at Butte. The problem thereby presented has never yet been clearly 
formulated, let alone any progress made toward solving it. 

Butte has long called itself the richest hill on Earth and this claim is re- 
garded by the geologists of the Anaconda Company as more than ever justi- 
fied (Lindforth, 1952). Three attempts have been made to account for the 
extreme localization of copper and zinc in that part of the bathylith in which 
Butte is situated. One explanation, by Weed, was that the granitic rock at 
Butte is richer in ferromagnesian minerals, a feature that favored precipitation 
of the copper. But granodiorite of the composition of the rock at Butte is 
common throughout the bathylith as a whole. Probably for this reason Weed 
omitted this explanation in his final report on the Butte district (1912). An- 
other explanation, much favored at Butte, is that the granodiorite at Butte 
was invaded by a porphyry, variously called quartz porphyry (Weed, 1897). 
Modoc porphyry, rhyolite porphyry (Weed, 1912, p. 40), or granite porphyry, 
which was directly responsible for the deposition of the ores. However, in- 
creasing knowledge of the areal geology and petrology of the bathylith shows 
that similar porphyry dikes occur in many places, but none of these places 
have copper deposits. A third explanation is that Butte is on an upward 
protuberance on the back of the Boulder bathylith, in short, a cupola, and 
this explanation would indeed be attractive if there were any field facts to 
support it. These attempts to explain the most remarkable economic feature 
of the bathylith emphasize the need of completing the areal mapping and 
petrologic study of the bathylith. 
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In conclusion, the present investigation, concerning itself mainly with 
the northern portion of the bathylith, has shown that the bathylith is not the 
result of the intrusion of a single magma in one pulse but was built up by 
the successive intrusion of at least five, probably six, bodies of magma. The 
problem of the emplacement of the bathylith is therefore not as Grout and 
Balk envisaged it in 1934 as the mechanics of intrusion of a single body of 
magma, but is the problem of how five different magmas in turn made room 
for themselves in the higher levels of the crust and built up a composite 
bathylith. 

These magmas after arriving in place were capable of developing facies, 
such as the granogabbro facies of the Unionville granodiorite. Furthermore, 
they have been altered by biotitization, whereby strikingly poikilitic dioritoid 
and gabbroid rocks were formed. In places also changes were induced by 
dissolving limestone. 

The successive major intrusive masses that make up the Boulder bathylith 
were emplaced in the normal order of increasing silica content. Whatever the 
process was that developed the successive magmas, it produced no change in 
the value of the alkalicity indices of the resultant rocks. All the rocks are 
normal lime-alkalic rocks. 

In those masses that do not fit into the normal scheme, such as the 
Scratchgravel Hills syenodiorite, the Priests Pass leucomonzonite, and the 
nepheline shonkinite, the alkalicity index falls below 55 and attains a low of 
12 for the nepheline shonkinite. 

The field evidence in regard to the nepheline shonkinite as well as the 
chemical evidence (the shonkinite contains 13 percent of lime) strongly sug- 
gests that the alkalic nature of the shonkinite was determined by the assimila- 
tion of limestone. Consequently, in the absence of field evidence for the origin 
of the syenodiorite and the leucomonzonite we must conclude that their alkalic 
character like that of the nepheline shonkinite is the result of the assimilation 
of limestone, but that the assimilation took place at greater depth, below the 
present level of exposure. 

As I mentioned in the beginning, it was the problem of the mode of em- 
placement of the Boulder bathylith that caused me to return to Montana. The 
bathylith and its enveloping rocks proved to be much more complex than was 
anticipated and detailed mapping therefore became necessary, and this map- 
ping led to the demonstration of the composite nature of the bathylith and 
unexpectedly to the discovery of striking field evidence in support of Daly’s 
theory of the origin of alkalic rocks. Many other interesting features were 
found such as the notable contact metamorphism, but these will require de- 
scription at a future time. 
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SEDIMENTARY CHARACTERISTICS OF DUST STORMS: 
Il. VISIBILITY AND DUST CONCENTRATION 


W. S. CHEPIL and N. P. WOODRUFF 


ABSTRACT. Analysis of some dust storms in Kansas and Colorado during 1954 and 
1955 indicates a relationship between visibility and atmospheric dust concentration when 
rules of Houghton are followed, Visibility varies inversely as some power of concentration, 
and concentration varies inversely as a certain power of height. The quantity of soil re- 
moved from any region for any storm or period of time can be estimated, 


INTRODUCTION AND ACKNOWLEDGMENT 

Visibility is often an indication of the concentration of dust in the 
atmosphere. Many weather observers in arid and semiarid regions record 
visibility to indicate the relative severity of dust storms. Such records are 
available for many years at different stations. Little information was found 
in the literature to indicate relationships between visibility and the actual 
quantity of dust carried by wind. The only available information of this 
nature is that of Langham, Foster, and Daniel (1938) giving visibilities as- 
sociated with various dust concentrations at 30 inches above the ground at 
Goodwell, Oklahoma, during 1936 and 1937. Analysis of these data is included 
in this report. Such data used in conjunction with available records of in- 
tensity-frequency of occurrence of windstorms should be of potential value in 
estimating losses of soil from wind-eroded regions. 

Estimations of soil losses also require information on the variation of 
dust concentration with height above the ground. Little of this type of in- 
formation is availab!e in the literature. While the data contained in this study 
do not include determinations to great heights, they do afford some oppor- 
tunity for at least speculating on the concentrations above the heights of 
measurement. The previous paper (Chepil, 1957) of this series indicated the 
proportion of dust lost from various soil classes by the process of deflation. 
This paper presents an analysis of concentration of dust in the atmosphere 
and of the approximate rates of its removal from wind-eroded areas. The study 
was conducted on some dust storms in western Kansas and eastern Colorado 
in 1954 and 1955. 

This paper is contribution no. 548, Department of Agronomy, Kansas 
Agricultural Experiment Station, Manhattan, and Agricultural Research Ser- 
vice, U. S. Department of Agriculture, Cooperative Investigations on the 


Mechanics of Wind Erosion. 


PROCEDURE 

Two grams of fine glass wool of 5-micron thickness of filaments packed 
in 1-1/8 inch inside diameter round aluminum tube was used to filter dust 
from the air. The apparatus consisted essentially of the glass wool filter, a 1- 
incu diameter connecting hose, an air meter and barometer to measure the 
volume of air intake as under standard temperature and pressure, an electric 
vacuum cleaner motor and fan to supply the necessary suction, a 110-volt 
generator, and a gasoline motor. Four filtering tubes were connected each to 
a separate manometer tube, and two such filtering tubes were connected to 
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each of the two vacuum cleaner units. Air intake in each tube was controlled 
by separate gas valves. The velocity of intake was the same as the velocity of 
the wind at corresponding height. A 1oultiple alcohol manometer and Pitot 
tubes were used to measure wind velocity at each location of the filter tube. 
The basic unit with the exception of the filter tubes, air meters, and barometer, 
was similar to that of Zingg (1951) for measuring total quantities of soil 
eroded by wind. The controlled intake device was used to gauge air volume 
of samples in 1954. Air meters were added for the 1955 studies. With the 
exception of filtering tubes the apparatus was also similar to that of Langham 
and others (1938) for measuring total weight of dust in the air. 

The assembly was mounted on a truck and hooded to give it and the 
operators partial protection from the dust storms (fig. 1). The filtering tubes 
were clamped to a vertical pole on the rear of the truck at various heights 
above the ground and facing into the direction of the wind. In 1954 the 
heights of measurement were 4, 6, and 8 feet, but in 1955 the heights of 
measurement were changed to 2, 5. 11. and 20 feet. 

The filtering tubes were dried in an oven at 110°C before and after each 
exposure to determine the weight of dust caught. The tubes were stoppered 
when not in use. The efficiency of the filtering tube for catching dust was com- 
pared with that of the impinger tube of Langham and others (1938). On the 
average the filtering tube caught 97.5 percent of the dust collected by the 
impinger. An appropriate quantitative correction was, therefore, made. The 
impinger tube has one serious disadvantage over the filtering tube in that it 
is impossible to determine from the sample of dust caught in the water the 
actual size of particles carried through the air. The filtering tube, on the other 
hand, merely trapped the dry dust which was shaken out readily and its 
equivalent size distribution determined. 


Fig. 1. Atmospheric dust-catching equipment as used in the field. Visibility in this 
case was 0.55 mile near Menno, Kansas, March 23, 1955. 
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The equivalent size distribution of dust particles was determined by sedi- 
mentation in carbon tetrachloride. Because this liquid is non-polar, it causes 
little or no breakdown of particles immersed in it. A weighed sample of dry 
dust first was wetted slowly by the liquid and then immersed in a definite 
volume of the liquid. The suspension was mixed thoroughly and allowed to 
settle for an appropriate period, depending on the equivalent diameter of the 
particles to be determined, depth of sampling, and temperature of the liquid. 
A 25 ce aliquot was then pipetted, using a uniform and appropriate air suc- 
tion from a pump. The aliquot was evaporated to dryness, and the residue 
weighed and expressed in percentage of the total weight of sample. 

Rules of Houghton (1945) for estimating daytime visibility were fol- 
lowed. Especially the following conditions were adhered to: 

1. Measurements were made only on cloudless days between 9 A.M. and 
P.M. 

2. Objects for visibility marks were as dark as possible, silhouetted 
against the horizon sky. Glittering objects were avoided. 

3. The sun was preferably not behind the observer but was in the field 
of vision. 

Dull-colored approaching vehicles coming within the field of vision often 
served as a measure of visibility. 


VARIATION OF DUST CONCENTRATION WITH HEIGHT ABOVE GROUND 

Relationships between height above the ground surface and measured 
dust concentrations for 1954 and 1955 and the 1954-1955 average are shown 
in figure 2. The 1954 curve is based on 10 measurements at heights of 4, 6, 
and 8 feet and the 1955 curve is based on 12 measurements at heights of 2, 
5, 11, and 20 feet. These relationships plot as a straight line on log-log paper 
and therefore may be expressed by a power equation of the generalized form 
a 


(1) 


where C is concentration of dust in milligrams per cubic foot at height y 
expressed in feet. The average value for constants a and y was found to be 
12.4 and 0.28, respectively (fig. 2). The constants varied little from one year 
to the other. It will be noted that the curves have been extrapolated to a 
height of 5282 feet. Any extension of these equations to heights greater than 
10 to 15 percent is extremely hazardous without some basis for doing so. In 
this connection a review of literature indicated that Schmidt as reported by 
Vanoni (1946) used an equation of the basic form 


y 
dy 
log -W f (2) 
a 


to express the concentration of dust in the atmosphere. If it is assumed that 
Es, the sediment transfer coefficient, is equal to E,,, the momentum transfer 
coefficient, Eg can be expressed in terms of the shearing force r, the depth of 
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flow ym, and any given height y, thus permitting integration of equation (2) 
to give 


Z 


Ym — a (3) 
where 
C = concentration at any height y 
C, = measured concentration at reference height a 
Ym = depth of flow (taken as one mile) 
T 
k 
where W = settling velocity of particle 
k = 0.4 
t = wind shear 
p = mass density of air 


In analyzing the present data, attempts were made to use this equation 
as a basis of extension of data. Fortunately it also plots as a straight line on 
log-log paper and once determined can be reduced to the basic form of equa- 
tion (1). Unfortunately, however, the degree of fit obtainable with this 
equation depends upon the determination of the mean equivalent diameter 
and the settling velocity of a representative sample of the suspended soil ma- 
terials and on an accurate determination of wind shear and depth of flow. 

In this study a rather thorough analysis of the equivalent diameter of 
dust particles in relation to height above ground surface was made for the 
data obtained during the 1955 season. This relationship and the equation 
expressing equivalent diameter d as a function of height y is shown in figure 
2. Since the mean equivalent diameter of the particles varies with height and 
also since a larger concentration of particles would be located nearer the 
ground, both a weighted mean equivalent diameter and a mean equivalent 
diameter were determined by dividing the total height of 1 mile into 11 
increments and integrating the two expressions for concentration and di- 
ameter, i.e., 


10.5 

(4a) 
and 

0.07 

(4b) 
thus giving 
Xd 
mean equivalent diameter = ———— = 0.0384 mm 
and 
XC 

mean weighted equivalent diameter = regent 0.0467 mm 
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Fig. 2. Measured and projected average equivalent diameter and concentration of 
dust at various heights for some 1954 and 1955 dust storms in Kansas and Colorado. 


Settling velocities W for the mean and the mean weighted equivalent diameters 
were taken as 12 and 15 cm/sec, respectively, from curves prepared by Rouse 
(1937) for the settling velocity of quartz spheres in air and water. 

The degree of fit of the Schmidt equation was found to vary with the 
diameter used. When the mean weighted equivalent diameter (0.0467 mm) 
was used, the equation indicated a greater decrease in concentration with 
height than the plotted data and it gave the closest fit to the 1955 data. Some 
further investigation revealed that when the mean equivalent diameter 
(0.038 mm) was used, the Schmidt equation was identical to the 1955 derived 
equation and gave a reasonable fit to the combined 1954-1955 data. It would 
not, however, fit the 1954 data unless a mean equivalent diameter of 0.025 mm 
was used. 

It is evident that there are some differences between the Schmidt equa- 
tion and the empirical equations derived in this study. However, the agree- 
ment that can be obtained with proper selection of mean equivalent diameter 
of particles, a privilege which is augmented by a general lack of definition of 
Schmidt's “mean diameter of representative sample” and the fact that one of 
the diameters giving a reasonable fit is a mean calculated from the data, would 
indicate that the empirical equations probably are reasonably representative 
of conditions to a one-mile height. The one-mile average depth of flow appears 


| 
| 


of Dust Storms: I1. Visibility and Dust Concentration 109 


to be reasonably valid. Occasionally dust clouds were reported to extend above 
12,000 feet. 

The total dust load in a cubic mile of the atmosphere can be found by 
determining the area under the height-concentration curves by integration. 
Values of the total dust load for each of the measurements made in 1954 and 
1955 are shown in table 1. The average total load for each year and for the 
combined years is also given. All of these values were computed on the basis 
of the measured concentration at a 6-foot height and the 1954-1955 average 
concentration-height relationship. In general, the total dust loads and rates of 
removal were less in 1955 than in 1954. This is not surprising since variations 
in both soil and climatic conditions during the two years would account for 
great differences in dust load and its rate of removal. 


RELATION BETWEEN VISIBILITY AND DUST CONCENTRATION 
Pertinent data pertaining to this phase of the study are given in table 1. 
Functional relationships between measured dust concentration at the 6-foot 
height (lower abscissa) and visibility are shown in figure 3. The equation of 
the curve derived by the method of least squares is: 


(5) 


where V is visibility in miles and C, is dust concentration in mg/cu ft at the 
6-foot elevation, Equation (5) is derived from data obtained in the present 
study and also from the data of Langham and others (1938) presented. in 
figure 3 primarily to show the general agreement between the two studies. 

Considerable data on visibility during dust storms are available from the 
Weather Bureau Stations. These visibility data probably could be used to 
determine dust concentrations. Information required would be a relationship 
similar to that of figure 3. The concentrations should be expressed, however, 
in terms of the quantity of material contained in a given volume of the at- 
mosphere. The upper abscissa of figure 3 converts the average curve to units 
of total dust load in tons per cubic mile. This conversion was made on the 
basis of measured concentration at 6-foot height and the average proportion- 
ality of concentration to height, C &% y~°-*8 
this relationship is: 


. The equation expressing 


29.5 


Co (6) 


where C,, is the concentration in tons per cubic mile of lower atmosphere and 
V is the visibility in miles. This equation provides a means of estimating dust 
concentrations from information on the visibility associated with a given dust 
storm. 

Visibility and dust concentrations are generally proportional to wind 
velocity, but there are major exceptions. For example, the highest wind 
velocity of 32.5 mph at 8-foot height at Cheyenne Wells, Colorado, on April 
27, 1955 actually was associated with a lower concentration than the wind of 
25.6 mph at Syracuse, Kansas, on March 10, 1954. Such results are to be 
expected since it is known that high wind velocity is not always associated 
with dust storms. The soil texture, the condition of the soil, the wind direc- 
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Fig. 3. Relation between visibility and dust concentration at 6-foot height. 


tion, the proximity to or position within the eroded area, and the previous 
intensity of the wind and erosional history of the region are some of the other 
factors that influence the concentration of dust in the air. 


SOIL REMOVAL IN RELATION TO DUST CONCENTRATION AND VISIBILITY 

The dust concentration-visibility relationships when combined with wind 
velocity data can be used in estimating the rate of soil removal from wind- 
eroded areas. The rate of removal past a given vertical area normal to wind 
direction would be equal to the concentration multiplied by the wind velocity. 
Since both the concentration and the wind velocity vary with height, the 
actual rate of removal for a given storm would be expressed as 


R=  X(C,u,)A (7) 
where R = rate of removal 
C , = increment concentration 
u, = increment wind velocity 
A = vertical cross-sectional area 


The most accurate method of summing up the increment velocity-con- 
centration products is to integrate the mathematical expressions for the varia- 
tion of velocity and concentration with height. The type of relationship 
between concentration and height is given by equation (1). The wind velocity- 
height relationship can be expressed as an exponential type equation plotting 
as a straight line on semi-log paper. This relationship is characteristic of a 
gradient wind which during the daytime hours in a fully turbulent atmosphere 
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was shown by data of Sutcliffe (1936) to extend at least to a one-mile 
height. Parkinson (1936) furthermore asserted that the presence of dust 
storms over the Great Plains area of the United States is usually associated 
with a high degree of atmospheric turbulence. The authors of this paper fully 
recognize the fact that under certain conditions the exponential relationship 
between velocity and height does not apply to the whole mile height. The 
assumption that it does apply to one mile is made as a matter of convenient 
but not at all unreasonable expedient. Even if one assumes that the exponential 
wind velocity-height relationship extends only 2000 feet and that there is no 
change in velocity from that height upwards, the computed rate of dust move- 
ment would not be more than 5 percent smaller than the rate based on the 
exponential relationship extending to a full mile height. The difference is small 
because both wind velocity and dust concentration vary relatively little between 
2000 and 5280 feet. 

Using the exponential relationships between concentration, wind velocity, 
and height, equation (7) may be written as the integral 


ad ] 

R= aA h + Bh ‘| dh (8) 

J clm 

where a = constant of concentration-height relationship 
B constant of velocity-height relationship 
m slope of velocity-height relationship (tangent of the 
angle the curve makes against the velocity axis) 

n slope of concentration-height relationship (0.28) 
A vertical cross-sectional area 
h height above ground 


Assuming that the concentration-height and the wind velocity-height 
relationships represent conditions to a one-mile height, inserting proper con- 
version factors and choosing the mile to extend from 2 feet to 5282 feet, but 
ignoring the 2-foot limit because of its negligible influence on the total rate, 
equation (8) can be integrated and simplified to give 


2.33 
R = 204a [= t B | (9) 
m 


where R is the rate of removal in tons/hour/vertical square mile area normal 
to wind direction. 

Average rates of removal through a vertical square mile area for the 
individual measurements of this study are given in table 1. These values were 
calculated from equation (9). 

While equation (9) permits estimates of rates of removal where both 
concentration and velocity profiles are measured, it does not permit calcula- 
tions using variables that are more readily available from Weather Bureau 
records. Two of the variables usually available are the visibility and the wind 
velocity at a given height. If equation (6) of this study is used to provide an 
expression for concentration in terms of visibility, then equation (9) can be 
expressed as 


29.5 2.33 
R B | (10) 
V1.2! m 
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Solution of this equation would require values of the slope m and the 
constant B of the velocity-height profile. These values would, of course, vary 
somewhat with different wind storms. If, however, wind velocity were meas- 
ured at one height and if the average value of Zo, i.e., the intercept of the 
velocity-height semi-log relationship for a given type of land terrain were 
known, values of B and m could be determined. Analysis of the 22 profiles 
obtained in this study indicated that on the average the value of Zo was ap- 
proximately 0.07 feet, which is probably a representative value for the flat 
terrain over which these profiles were measured. Wind velocity at 60 feet was 
estimated from the velocity data up to 20 feet simply by extending the straight 
line curve of the plotted velocity-log height relationship to 60 feet. Values of 
B and m in terms of this value and the wind velocity at a height of 60 feet 
are as follows: 


Veo (mph) m B 
20 0.1477 7.82 
25 0.1182 9.77 
30 0.0984 11.74 
35 0.0844 13.68 
40 0.0738 15.65 
45 0.0656 17.61 
50 0.0591 19.54 
55 0.0537 21.51 


60 0.0492 23.47 


These values used in conjunction with visibility measurements substituted into 
equation (10) would give an estimate of the rate of soil removal. 

The next step in this type of analysis would be to determine the quantity 
of soil removal for any storm or period of time. While equation (10) pro- 
vides a means of estimating rates of removal, it cannot give a total or actual 
removal without some knowledge of the width of the eroded area normal to 
wind direction and duration and number of dust storms occurring in a given 
length of time. Considerable data are available from Weather Bureau records 
on visibility, wind velocity, and duration of dust storms. This paper merely 
indicates how these data might be used in estimating the quantities of soil 
material removed from wind-eroded regions. 
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EXPERIMENTAL DETERMINATION OF 
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* Paper No. 147 published under the auspices of the Committee on Experimental Ge- 
ology and Geophysics and the Division of Geological Sciences at Harvard University. 
ABSTRACT. The reaction, albite + nepheline — 2 jadeite, has been studied experi- 
mentally at pressures between 10,000 and 25,000 kg/cm’, and at temperatures betwee: 
600° and 1200°C. The pressure system is an enlarged version of Bridgman’s double press, 
in which a conical high-pressure chamber is given external support by heavy rings. The 
pressure medium is nitrogen, Temperatures are produced by an internal wire-wound 
furnace, in which are heated small gold or platinum capsules containing the charges, 
Quenching is obtained when the furnace power is cut off. Charges of natural jadeite, 
natural albite plus synthetic nepheline, natural analcite, and synthetic jadeite glass have 
been heated for periods up to several hours under pressure, quenched, and examined with 
the aid of an X-ray diffractometer and under the microscope. 


The equilibrium 1 line for the reaction is considered to be satisfactorily represented 
by the equation, P = 1000 + 18.5 T, where P is the pressure in kg/cm*, T the tem- 
perature in degrees centigrade. The uncertainty of the intercept for T = 0 is + 1000 
kg/cm*, the uncertainty of the slope, + 1 (kg/cm*)/degree. The experimental line is 
in reasonable agreement with the thermochemical data. Jadeite crystallizes on the high- 
pressure side of this line from glass of jadeite composition, or from an albite-nepheline 
mixture. The position of the jadeite melting curve has been roughly determined. 


The discussion relates these results to the question of the high-pressure breakdown 
of the feldspars, and to the possibility of accounting for the Mohorovicic discontinuity as 
a boundary for such transformations, 


INTRODUCTION 


The conditions of formation of the mineral jadeite have attracted much 
attention in recent years. Petrologically, jadeite is interesting not so much 
for the relatively rare occurrence of the pure mineral as for its presence as 
a component in the pyroxene, omphacite, which is found in eclogites. Eclogites 
have been thought to have formed at the highest pressures attained during 
regional metamorphism (Eskola, 1939), and until recently, despite many at- 
tempts, jadeite had not been synthesized. The highest pressure employed in 
the unsuccessful experiments, which have been reviewed by Yoder (1950b), 
appears to have been no greater than 4000 bars. Within the last few years 
jadeite has been produced experimentally at pressures greater than 10,000 
bars in at least three laboratories, and it has become possible to begin a syste- 
matic exploration of its field of stability. 

Jadeite (NaAlSi,O,) has a composition intermediate between that of 
albite (NaAISi,Os) and of nepheline (NaAISiO,.) But whereas the latter 
minerals have framework structures and densities of about 2.6, jadeite has 
a pyroxene structure and a density of about 3.3; its compressibility is com- 
parable with that of olivine (Adams and Gibson, 1929; Yoder and Weir, 
1951.) The “jadeite problem” may thus be considered as part of a problem 
of much wider significance: the feldspar problem. The feldspars are with 
little doubt the principal minerals of the Earth’s crust, that is, of the outer- 
most shell which is bounded internally by the Mohoroviéié discontinuity of 
the seismologists. Below this discontinuity the seismic velocities assume values 
which are definitely incompatible with the velocities in rocks containing 
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major proportions of any of the feldspars. There are *wo conceivable explana- 
tions: (1) a discontinuous change of composition, from gabbroic to ultra- 
mafic, with a drastic impoverishment of the subcrustal material with respect 
to alkalies and alumina; or (2) a change of crystalline phase, from feldspar 
to jadeitic pyroxene and other high-density minerals, especially garnets, with 
little immediate change of composition. Though the first possibility has won 
wide acceptance, the second, which is essentially a form of the “eclogite 
hypothesis,” still deserves careful attention. This hypothesis originated with 
Fermor (1912, 1913, 1914, 1938) and has been advocated by Joly (1930, 
p. 53-54) and by Holmes (1945, p. 507.) The following study of the jadeite 
equilibrium is, of course, merely a first step in the experimental investigation 
of the eclogite equilibrium; further work must include at least lime and 
magnesia as additional components. 

The thermochemical data relating to the jadeite equilibrium have been 
studied in several recent papers. The heats of formation of the required com- 
ponents have been determined by Kracek, Neuvonen, and Burley (1951), 
and the entropies by Kelley, Todd, Orr, King, and Bonnickson (1953.) These 
writers, as well as Yoder and Weir (1951), and Adams (1953), conclude that 
jadeite is stable with respect to albite plus nepheline at ordinary pressure and 
temperature; the equilibrium curve for the reaction ab + ne = 2jd has been 
computed by Adams and by Kelley et al. 

The present work was begun in late 1952, and jadeite was first identified 
with certainty among the synthetic products in May, 1953. Unknown to us, 
a program of high-pressure synthesis had been in progress for several years 
at the Norton Company of Worcester, Masssachusetts, and jadeite, as well as 
pyrope, kyanite, and a new dense polymorph of silica (Coes, 1953, 1955). 
had been produced at pressures between 20,000 and 45,000 bars; the synthesis 
of jadeite was first accomplished in September, 1948.’ The breakdown of 
analcite to jadeite, and to albite plus nepheline, at high pressures has been 
reported by Griggs, Kennedy, and Fyfe (1955a, 1955b). Kyanite and pyrope 
have also been crystallized in the course of the present work (Robertson, 
Birch, and MacDonald, 1955); their equilibrium relations will be discussed 
in other publications. 
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EXPERIMENTAL DETAILS 
Pressure system.—tThe pressure system is closely patterned after a smaller 
one designed by Bridgman (1938) for use with liquids; the linear dimensions 
of our equipment are approximately 50 percent greater than those of Bridg- 
man’s apparatus. Thus the bore of the high-pressure chamber is 3/4 inch 
instead of 1/2 inch; this enlargement provides barely adequate space for an 


Fig. 1. Arrangement and principal components of pressure system: A, tie rods and 


plates of press; B, lower ram; C, upper ram; D, supporting rings; E, high-pressure 
chamber. 
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internal furnace. The requirement of temperatures of the order of 1000°C 
excludes the use of liquids, and the pressure medium was either argon or 
nitrogen. 

The principal components of the high-pressure system are (fig. 1) a 
framework (A) of plates and tie rods, a lower hydraulic ram (B), an upper 
hydraulic ram (C), the tapered high-pressure cylinder (E) and its support 
rings (D) (see also fig. 2.) Two pumps capable of producing liquid pressures 
of about 20,000 lbs/in? supply the two hydraulic systems. The upper ram 
forces a tungsten carbide piston into the bore of cylinder E to produce the 
internal pressure; at the same time, the lower ram forces E into the support 
rings D, to produce an external pressure over the tapered surface. The ratio 
of internal pressure to external pressure was usually maintained at about 3/1. 
Before pressure can be built up with this system, a preliminary compression 
of the gas raises the pressure from that of an ordinary commercial cylinder to 
about 20,000 lIbs/in?. 

A number of details concerning packings, clearances, materials, and so 
on, were varied from time to time in efforts to eliminate leaks, and to improve 
the life of the parts and convenience of access. These have little general interest 
and will not be described here but they will be recorded in a separate report 
to the Office of Naval Research. It is perhaps worth recording, however, the 
most common type of mechanical failure. This was a rupture of the high- 
pressure cylinder on a surface approximately normal to its axis and several 
inches from its small end, often accompanied by a radial fracture down the 
length. This type of rupture. which sometimes occurred at pressures as low 
as 12,000 bars, appears to be related to the discontinuity of internal pressure 
at the piston packing; below this packing, the pressure has its high value; 
above the packing, it is zero. At the same time, the external surface of the 
tapered cylinder is acted upon by the supporting pressure, which is not uni- 
form but reaches its highest value near the small end. Improvement resulted 
from reducing the external pressure in the uppermost support ring by giving 
this ring a small amount of initial clearance, but it cannot be said that the 
difficulty has been finally eliminated. 

The furnace (F in fig. 2) is assembled in a thin-walled steel tube which 
slips into the bore of the high-pressure chamber. The bore of the porcelain 
furnace tube is 1/4 inch, the length of winding about 3 inches on a diameter 
of 5/16 inch. The space between the winding and the steel tube is filled by 
a sleeve of “Lava” (mostly pyrophyllite, marketed by the American Lava 
Corp., Chattanooga) which serves as electrical and thermal insulation. Within 
the furnace bore is a multibore porcelain tube carrying thermocouples. The 
upper part of the furnace is closed by a solid plug, and all interior empty 
spaces filled with powdered alumina, in order to reduce convection. The 
furnace winding is connected at one end to the steel sleeve, which is grounded; 
at the other, along with the thermocouples, to insulated, pressure-tight stems 
in the bottom closure. In a convenient variation, a miniature socket connector 
is used as part of the furnace tube and makes an easily disconnected junction 
with the leads of the bottom closure. Also at the base of the furnace tube, 
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D 
Cc 
A 
\ 
Fig. 2. Internal arrangements of high-pressure system: A, high-pressure cylinder; 


C, supporting rings; B, lower closure, with electrical leads; D, piston mushroom, with 
packing; E, piston; F, furnace. 


below the heated region, is the manganin pressure gauge, and in some ar- 
rangements, a copper resistance thermometer. 

The furnace winding is of nichrome, kanthal, or platinum-rhodium wire 
or strip, usually designed to have a hot resistance of about 13 ohms. The 
power input is controlled by a variable transformer supplied by a constant 
voltage transformer with a capacity of 1 KVA. 

Measurement of pressure and temperature-—Though rough measurements 
of pressure may always be made from the ratio of areas of piston and ram 
and the pressure on the ram, the uncertainties resulting from friction and dis- 
tortions may be large; the pressures given in this work depend upon the 
change of electrical resistance of a manganin coil exposed to gas pressure. 
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The high-pressure scale depends upon fixed points established by Bridgman 
(1940), of which we have used the freezing pressure of mercury, 7640 kg/cm? 
at O°C and 12,000 kg/cm? at 22°C, and the transition pressure for the bismuth 
modification I to Il, 25,800 kg/cm? at 22°C.* These pressures are probably 
accurate to about one part per thousand. With the aid of these fixed points, 
Bridgman has shown that the resistance of manganin wire is not quite linear 
with pressure: if the coefficient of a manganin coil is determined at the freez- 
ing point of mercury and then used, by linear extrapolation, to measure the 
bismuth transition pressure, the error is about 2 percent. This would not be 
a significant error for the present work; it can, moreover, be eliminated by 
the use of a second-degree relation between resistance and pressure. The 
principal uncertainty of pressure measurement in this work arises from the 
effect of temperature upon the resistance of the manganin gauge coil. 

Between 20° and 30°C, the resistance of manganin passes through a flat 
maximum with respect to temperature, and in ordinary installations of man- 
ganin pressure gauges it is rarely necessary to correct for small changes of 
ambient temperature. In the present system, however, the gauge may reach 
temperatures as high as 100°C, which) depend upon the power input and the 
duration of heating. As the coil temperature rises at constant pressure, its 
resistance decreases, giving the appearance of a drop in pressure. It became 
necessary to determine the temperature of the coil, and also to examine the 
effect of temperature upon the resistance-pressure characteristic of manganin. 

The pressure gauges were 50-ohm coils of no. 40, double-nylon-covered 
manganin wire recently obtained from the Driver-Harris Company, wound 
non-inductively on Teflon spools, and protected by Teflon sleeves. The sensi- 
tivity of resistance measurement corresponded to a pressure difference of 
10 kg/cm*. A number of these coils were calibrated against the freezing point 
of mercury at 0°C at various times; the pressure coefficient is 2.31-10~° per 
kg/cm®* at room temperature, with variation of about 0.5 percent. (Coils made 
of enameled manganin wire obtained from the same manufacturer about 
twenty years ago show a coefficient of about 2.5-10—° kg/cm*.) Several coils 
were then heated to temperatures up to about 125°C at 1 atmosphere, and 
also at 7500 kg/cm*, measurements of resistance being taken as function of 
both temperature and pressure. The results are summarized in table 1. The 
effect of temperature is slightly smaller at 7500 kg/cm? than at 1 atmosphere, 
with the corollary that the pressure coefficient is slightly larger at 100°C than 

TABLE 1 
Effect of Temperature on the Electrical Resistance 
and Pressure Coefficient of Resistance of Manganin 


re R/R:z Increase of pressure coefficient, percent 
25 1.0000 0 
50 0.9998 - 
75 0.9991 0.8: 
100 0.9979 1.6; 
125 0.9965 2.7; 


* The unit of pressure for the experimental work is the kg/cm*, which equals 0.98 bar; 
the difference between these two units is of no consequence for the present work. 
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at 25°C. The effect of changing the temperature, at constant pressure, from 
25°C to 100°C is equivalent to a pressure change of approximately 900 
kg/cm*. An error in the temperature of the coil of 5°C at 100°C is equivalent 
to a pressure change of roughly 100 kg/cm*. The cumulative errors may 
amount to as much as several hundred kg/cm? at the highest pressures. 

For comparison with table 1, Bridgman (1938, p. 170) gives an increase 
of coefficient of 0.6 percent between 30°C and 75°C; an earlier and probably 
less satisfactory determination (Bridgman, 1935, p. 79) gives an increase of 
1.3 percent between 0 and 50°C, 2.0 percent between 0 and 95°C. There is 
some evidence that the temperature coefficient of the pressure coefficient is 
more variable than the pressure coefficient itself. 

The temperature of the pressure gauge in our high-pressure system was 
obtained from the change of resistance of a copper coil mounted close to the 
pressure coil and having the same external dimensions; it was supposed, for 
this correction, that the slope of the resistance-temperature curve of copper 
is unaffected by pressure, in accordance with Bridgman’s determination 
{1931, p. 262; 1938, p. 174) of the small effect of temperature upon pressure 
coefficient. A rough value of the coil temperature may also be obtained from 
the temperature of the base of the tapered cylinder. 

Temperatures in the furnace were determined by thermocouples, usually 
with at least two junctions close to the charges. The temperature is far from 
uniform; it varies from as much as 1200°C or so at the hottest region to about 
100°C at the cylinder wall, less than 1 cm distant radially, and to about the 
same temperature at the ends of the furnace, 3 or 4 cm distant axially. The 
best location for the charges was determined by trial. Most of the uncertainty 
in the temperature of the charges arises from the sharp variation of tempera- 
ture along the furnace; unless the charges were accurately centered upon the 
locus of maximum temperature, the variation along a single charge might 
reach 50° or more. By placing the charges between two thermocouple junc- 
tions, we attempted to bracket the actual temperature, but the uncertainty 
remains appreciable. 

Chrome-alumel couples were first used, with leads of chromel and alumel 
through the packing of the bottom closure to an external cold junction. The 
effect of pressure on the thermal electromotive force of this couple has been 
shown to be negligible up to 4000 kg/cm* and 600°C (Birch, 1939); the 
pressure has been extended to 7000 kg/cm? in unpublished work. It has been 
assumed that the pressure effect is still relatively unimportant at the higher 
pressures of the present study. Under ordinary conditions, the chromel-alumel 
couple may be used to some 1100-1200°C, but there is evidence of some kind 
of deterioration in high-pressure nitrogen at temperatures of 900°C and 
higher. This became noticeable as an apparent decrease of temperature at 
constant power input and constant pressure, and it has been tentatively 
ascribed to nitriding of the alumel wire. The effect disappeared when plati- 
num-platinum-10% rhodium couples were substituted for chromel-alumel. A 
pressure effect has been noted for this couple, amounting to an apparent de- 
crease of about 12°C for 20,000 kg/cm? and 1000°C on linear extrapolation 


from the measurements at lower temperatures and pressures (Birch, 1939). 
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A further uncertainty was introduced by the use, with this couple, of copper 
leads through the bottom closure. This made the connection from couple to 
copper inside the pressure chamber the effective cold junction, for which the 
temperature was taken as equal to that of the copper coil described above. 
This particular correction has been eliminated in later work by the use of 
platinum and platinum-rhodium leads through the packing. The major un- 
certainty as to temperature arises from non-uniformity of the temperature 
distribution, as mentioned above. 

Materials—The determination of the jadeite equilibrium line depends 
chiefly upon the crystallization of synthetic glass of the composition of jadeite, 
prepared by Dr. J. F. Schairer of the Geophysical Laboratory of Washington. 
Several other materials used for various purposes were as follows: 

Nepheline, synthetic, also prepared by Dr. Schairer. An analysis of this 

material is given in Kracek and others (1951.) 

Albite, Amelia, Virginia. An analysis of another sample from this locality 
is given in Kracek and others (1951). 

Quartz, Lisbon, Maryland. This is described in Kracek and others 
(1951). 

Analcite, Cape Blomidon, Nova Scotia, A sample of very pure analcite 
from the Holden Collection, Harvard University. 

Gibbsite, synthetic, Norton Company, Worcester, Mass. This is a com- 
mercial product; it was heated for 3 hours at 110°C in order to remove 
excess water, before use. 

Jadeite, Burma, no. 10464 of the Holden Collection, Harvard University. 

Procedure.—Finely ground samples weighing about 50 milligrams were 
loaded into 5/8-inch lengths of gold or platinum tubing, 1/8-inch in diameter ; 
thes: capsules were either left unsealed, with crimped ends, or sealed by arc- 


we srding to procedures developed by Yoder at the Geophysical Lab- 
Orato Three or more of these capsules could be packed in the furnace, 


between thermocouple junctions. 

As in most studies of silicate equilibrium relations, advantage has been 
taken of the sluggishness of the reactions: after an exposure to the highest 
temperature of an hour or so (see table 2), the charges were quenched to 
ordinary temperatures while still under pressure by the simple method of 
cutting off the power to the furnace. The transfer of heat in the dense nitrogen 
is so rapid that the temperature falls to about 100°C (or the temperature of 
the steel cylinder) in a matter of seconds, and it is assumed that no significant 
change of phase takes place in this time. The pressure was then reduced, and 
the charges removed for examination. 

In the beginning of an experiment, the pressure was raised with the 
furnace cold, to a pressure somewhat below the desired final pressure. The 
temperature was then raised to the desired point over a period of about 1 
hour; with the increase of temperature of the nitrogen, the pressures also in- 
creased, usually by 1000 to 2000 kg/cm’. The pressure and temperature were 
then kept as constant as possible at the maximum values, which are the ones 
indicated in table 2. It is assumed that, whatever reactions occur during the 
heating period, the final products correspond to the conditions of maximum 
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temperature. This assumption appears to be justified by the results. When the 
products are assigned positions in the pressure-temperature diagram according 
to the highest pressures and temperatures, there is a clear separation of fields, 
which could hardly result if the rate of heating or duration of heating had a 
significant effect. It may be noted furthermore that the approach to the high 
temperature is in all cases made by way of the jadeite field, as determined 
below, and similarly, on quenching the charge is returned to the jadeite field. 
Charges crystallized in the abite-nepheline field, however, showed no jadeite 
in the product. 

Identification usually depended upon comparison of X-ray diffraction 
patterns taken with a Norelco high-angle recording diffractometer. In some 
cases this was supplemented by microscopic determination of refractive index 
with the aid of immersion oils. Because of the purity of the starting materials, 
the number of possible products was small, and it was usually possible to ac- 
count for all significant lines in the diffractometer pattern. The amount of 
crystallization, when the initial ingredient was glass, varied from zero to about 
50 percent, with individual grains up to about 0.1 mm but usually smaller 
than 0.05 mm. The whole charge was generally used for the X-ray sample; it 
was then scraped from the mount and saved for possible future examination. 


A sample of the diffractometer records is shown in figure 3, with a com- 


, 


co 
8 
30 


40 


on 


ad 
a 


30 


BURMESE JADEITE JADEITE + CORUNDUM, FROM 


JADEITE GLASS + GIBBSITE 
20 


Fig. 3. Diffractometer patterns of Burmese jadeite and of synthetic jadeite and 
corundum formed from jadeite glass and gibbsite at about 800°C and 20,000 bars. The 
ordinate is 29 (CuKa), The corundum lines are labeled CO. 
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parison between the patterns of natural Burmese jadeite and of synthetic 
jadeite crystallized from synthetic jadeite glass plus a small amount of gibbsite. 


RESULTS 

The fields of stability can be roughly outlined with the aid of observa- 
tions of the behavior of natural minerals. On heating in the (ab + ne) field, 
natural jadeite breaks down to these minerals; similarly, on heating in the 
jadeite field, a mixture of albite and nepheline begins to form jadeite. The 
reactions are slow, especially below about 900°C, and never reach completion 
in the present experiments, even with the aid of such catalysts as water, gibb- 
site, or sodium fluoride. The results are significant only if a change of phase 
can be definitely established; this requires some minimum amount of trans- 
formation, possibly 5 to 10 percent. In about one-half of the runs starting 
with the crystalline phases, in which a change might be expected, the trans- 
formation could be observed. Only these significant results are plotted in 
figure 4; they suffice to fix a low-pressure limit for the reaction, given by 
P = 19 T (P in kg/cm’, T in degrees Centigrade.) They are consistent with 
a line drawn through 0°C and zero pressure, having a slope of 20 bars/deg 
or 50 deg/kilobar, as given by the thermodynamic data. On the other hand, 
because of the wide separation of the significant points at the lower tempera- 
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Fig. 4. Transformations of natural minerals in the temperature-pressure plane of 
jadeite composition, The starting materials were either Burmese jadeite or Amelia albite 
plus synthetic nepheline, as indicated. For details, see table 2. 
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tures, the high-pressure limit is given by P = 8000 + 13 T, so far as these 
points alone can define it, 

A closer determination may be obtained from the runs in which the 
starting material was synthetic glass of jadeite composition. Crystallization of 
this glass takes place more readily than does recrystallization of the minerals, 
and X-ray patterns indicating crystalline phases were obtained in nearly all 
cases. The amount of crystallization is sensitive to catalytic action, however, 
the best results being obtained with unsealed capsules containing initially 
about 15 percent of gibbsite which breaks down at a moderately high tem- 
perature to corundum and water. Free water was introduced on a number of 
occasions but, unless the capsules were sealed, it was probably lost too soon 
to be of value; if sealed in, on the other hand, very small amounts of water 
were sufficient to produce a low-melting silicate-water melt or solution which, 
on quenching, yielded only glass. It is assumed that water does not enter into 
the composition of the crystalline phases and, consequently, that the position 
of the equilibrium curve is not affected by its presence. At temperatures be- 
low about 600°C, however, analcite may be formed, and water becomes a 
true component in the reactions, analcite = jadeite plus water, and analcite 
= albite plus nepheline plus water. 
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Fig. 5. Transformation of synthetic glass of jadeite composition. The line marked 
“theoretical” is the equilibrium curve for the reaction, ab + ne = 2 jd, according to 
the calculations of Kelley et al. (1953). For details, see table 2. 
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The results for the runs in which glass was the starting material are 
plotted in figure 5, according to the crystalline phases identified after quench- 
ing from the indicated pressures and temperatures. In addition to the cases 
in which only jadeite, or only albite and nepheline (and corundum, when 
gibbsite was used) could be distinguished, there are a number of instances 
where jadeite, albite, and nepheline were found together; with one or two 
exceptions, these were cases where no catalyst had been added, or where only 
water had been added and had probably been lost. Nearly all of the corre- 
sponding points lie in the jadeite field, as already established by the reaction 
of mineral albite and nepheline. It seems probable that the transformation of 
the glass to crystalline jadeite passes through a stage of metastable crystalline 
albite plus nepheline and that, if the reaction is not sufficiently rapid, the 
quenched material shows both phases. Of the runs with gibbsite, only one 
(no. 35) shows this effect. If this explanation is correct, then we may consider 
that all of the runs which show crystalline jadeite in the product lie in the 
jadeite field. The equilibrium line is then fixed within relatively narrow 
limits; the best graphical solution falls close to the line, P = 1000 + 18-5T. 
The intercept at T = O°C, found by linear extrapolation, is subject to an un- 
certainty of perhaps +1000 kg/cm*; the position is best determined in the 
neighborhood of 900°C to 1000°C, where the uncertainty may be several 
hundred bars, or about 20 degrees, and there are no useful points below about 
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Fig. 6. Transformation of natural analcite and synthesis of analcite. For details, 
see table 2. 
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600°C. The uncertainty of the slope is about +1 bar/deg. Presumably, in 
the absence of water, the portion of the equilibrium line below 600°C could 
be realized given sufficient time or suitable non-reacting catalysts. 

Additional support for this line is given by the experiments in which the 
starting material was natural analcite. The results are plotted separately in 
figure 6. Three regions may be distinguished: in the first, analcite breaks 
down to albite plus nepheline, with the albite lines usually much the stronger 
in the diffractometer pattern; in the second, jadeite appears along with albite 
and possibly nepheline; finally, in the third region, below 600°C, analcite is 
unchanged, while either jadeite plus water, jadeite glass plus water, or albite 
plus nepheline plus water, may be transformed to analcite. The first region is 
evidently consistent with the (ab + ne) field already found; the second is 
consistent with the jadeite field already found, but it appears that the break- 
down of analcite, like the crystallization of the glass, may pass through a 
metastable phase of ab + ne. In the jadeite field, this in turn begins to trans- 
form to jadeite, and thus these runs may also be used to narrow the region 
of uncertainty of the equilibrium line. The third field, in which analcite is 
stable, has been little studied in the present work, but has been outlined by 
Griggs, Kennedy, and Fyfe (1955a; 1955b.) The equilibrium between analcite 
and albite plus nepheline plus water has also been investigated by Yoder 
(1954.) It should be noted that in the present work the reported pressure is 
the total pressure; the partial pressure of water is not usually known. 

The observations are summarized in table 2. Charges showing no change 
have not been listed. Phases showing only relatively weak lines have been dis- 
tinguished by parentheses. The uncertainties assigned to the pressures indicate 
the range of variation of pressure during the exposure to high temperature; 
the absolute pressures are subject to the errors discussed above. The tempera- 
tures include the uncertainty arising from temperature gradients as well as 
fluctuations during the run but are also subject to additional errors of absolute 
value. Many of these runs are for conditions well removed from the boundary 
and contribute little to the determination of its position, except as part of the 
process of narrowing down the limits. Only those closest to the boundary have 
been plotted in figure 5. The point given by Coes (1955) as representing the 
“best conditions” for the formation of jadeite, at 900°C and 20,000 bars, is 
indicated by a cross; this point lies well within the jadeite field as determined 
by the present work, about 2000 bars above the equilibrium pressure for 
900°C. 

Dry charges of natural jadeite and of albite plus nepheline were melted 
on two occasions to give some preliminary information on the liquidus rela- 
tions at high pressure. In the first, the pressure was 21,350 kg/cm*; the tem- 
perature was first held at 1280°C, then lowered to 1190°C and held for 20 
minutes, then quenched. Albite and nepheline formed from the glass made 
from the natural jadeite; the ab + ne mixture was still all glass. Accordingly, 
the higher temperature has been plotted in the liquid field, the lower in the 
ab + ne field. In the second case, on quenching from 1350°C and 25,100 
ke/cm*, the capsules, filled as before, contained only glass; and hence this 
point has been plotted in the liquid field. Dry natural jadeite did not melt 
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TABLE 2 
Summary of Experimental Results 
Abbreviations 


ab—albite jd—jadeite 
ana—analcite jd gl—jadeite glass 
co—corundum ne—nepheline 
gi—gibbsite wa—water 


(S)—sealed capsule 


Time Pressure Temperature Crystalline 
Run hours kg/cm* Cc Charge products 
77 2.0 9,000 + 30 535+15 jd + wa (S) jd + (ana) 
75 1.4 9.350 + 30 610+10 ana ab + ne 
jid + wa ab + ne + ana 
72 2.0 10,300 = 100 500+ 10 ab + ne + wa ana + ab + ne 
71 1.6 10,700 + 50 630+ 10 ana (5S) ab + ne 
80 1.5 12,000 + 30 575+ 35 id gl + wa ana 
20 2.0 12,200 + 200 775 +25 jd ¢l + NaF (ab) + ne 
19 1.5 12,400 + 100 1015+ 30 id gl ab + ne 
96 3.5 12,700 + 200 610 +20 id gl + gi ab + ne 
53 1.3 12.800 + 100 1050+ 25 id jd + ab + ne 
54 08 12,800+100 1025+20 id jd + ab + ne 
id gl ab + ne . 
74 2.0 13,400 +50 §25+15 ana jd + ab 
21 1.0 13,300 + 300 675 +25 ana ab + ne 
jd gl + gi (S) ab + ne 
24 1.4 13,000 + 300 775 +25 jd gl + NaF (ab) + ne 
jd + gi (jd) + ab + ne + (co) 
ana ab 
id gl + gi ab + ne + (co) 
70 1.7 14,000 + 50 610+ 10 ana jd + ab + (ne?) 
25 3.0 13,900 + 500 750+50 id gl + gi ab + (co) + (gi) 
jd + gi jd + (ab) 
ana jd + ab + (ne?) 
63 1.5 15,400 + 100 825 +25 jd gl + wa ab + ne 
id + wa jd + ab + (ne) 
ana ab + ne 
67 1.6 15,700 + 100 615+ 10 ana jd + ab + (ne) 
ana (S) jd + ab + (ne) 
600 + 10 jd gl + wa jd + (ab + ne) 
ab + ne + wa jd + ab + ne 
60 2.0 15,850 + 50 825+ 10 ana ab + (ne) 
61 2.0 16,100 + 100 825+ 10 ana ab + ne + jd 
jd gl + wa ab + ne + jd 
69 0.5 16,200 + 200 1150 +20 jd gl ab + ne 
id jd + ab + ne 
59 23 16,500 + 100 815+15 ana jd + ab + (ne) ; 
32 15 16,700 + 100 950 + 20 jd gl + gi (S) ab + ne + co 
45 1.8 16,700 + 200 875+15 jdgl (S) ne + (a5) 


58 15 18,000 + 300 920 + 20 jd gl ab + ne 


Time 
Run hours 


68 


2.5 


w 


0.1 
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Pressure 
kg/cm? 


18,300 + 200 
18,300 + 200 
19,000 + 50 


19,200 = 200 
19,400 + 50 


19,450 + 50 
19,500 + 50 
19,000 + 500 
19,400 = 300 
19,200 + 700 
19,300 + 600 


19,600 + 400 
20,100 + 100 


21,250+ 50 


21,350 + 350 


21,350 + 50 
21,800 + 200 
21,200 + 800 


23,900 + 30 


25,500 = 100 


26,000 + 50 


TABLE 2 (Continued) 


Temperature 


950 + 20 


1035 +15 
1080 + 20 


1020+ 10 
1055 + 25 


1190 +25 
995 +10 
900 + 20 


1080 + 50 


1175 +50 


1085 + 30 


129 


Crystalline 
Charge products 
jdgl + gi jd + co 
jd gl + NaF jd 
jd gl + gi (S) jd + co 
ana (S) jd + ab + (ne) 
jd gl + gi (S) jd + co + ab 
ana (S) ab + (jd + ne) 
ana (S) jd + ab + (ne) 
jd gl jd 
ana jd + ab + (ne?) 
jd gl + gi (S) jd + co 
ana (S) jd + ab + (ne) 
jd gl + gi (S) jd + co 
jd gl + wa (S) (ab) + ne 
jd gl + gi (S) (ab + co) 
jd gl ab + ne 
idgl (S) ab + ne 
jdgl + wa (ab) 
jd gl + gi (S) jd + (co) 
ana (S) jd + ab 
jdgl + wa ab + ne 
jd gl ab + ne 
jdgl + gi (S) jd + co 
jdgl + wa ab + ne 
jd + wa jd + ab + (ne) 
jdgl (S) ab + ne + (jd) 
jd gl jd + ab + ne 
ab + ne jd + ab + ne 
jd + wa ab + (ne) 
jd gl jd + ab + ne 
jdgl+gi (S) jd 
ana (S) jd + ab + (ne) 
jdgl + wa jd + ab + ne 
id gl jd + ab + ne 
ana jd + (ab) 
ab + ne ab + ne +jd 
jd gl jd + (ab) 
jd jd (not melted) 
jd gl jd + ab + ne 
ab + ne + (jd) 


ab + ne 


in 1.4 hours at 1175° and 25,500 kg/cm?. The very tentative construction of 
the liquidus and solidus curves in figure 5 is based on these results, plus those 


for one atmosphere given by Greig and Barth (1938.) 
The melting relations with water at high pressures remain to be worked 


out. The results with sealed capsules indicate reductions of melting point of 
several hundred degrees, but the water content has not been determined. A 


34 O15 840+ 60 
3 OS 950 + 20 
36 (15 950 + 20 
SJ 20 935 + 30 
28 «10 920+ 20 
38S 975 +25 
0 1.0 1000 + 100 
4 15 1000 + 25 
304 900 + 20 
3 1000 + 50 
37 
56 1.3 
81 0.3 
31 
6 14 
= 1.4 
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beginning has been reported by MacKenzie (1954) for water pressures up 
to 1000 bars. 

In many of the runs listed in table 2, capsules were included with charges 
having compositions of other systems, such as kyanite and pyrope. Both of 
these minerals have been crystallized, but further work is required for the 
delineation of their fields of stability. 


DISCUSSION 


Thermochemical considerations——The thermodynamic properties of 
jadeite, nepheline, and albite are, for silicates, exceptionally well known. A 
comparison of the equilibrium relations derived from thermochemical data 
with the experimental results is useful in several ways. First, the theoretical 
relations refer to equilibrium conditions, while the experimental curve is 
derived from short-time experiments in which equilibrium may not have been 
reached. The equilibrium curve calculated by Kelley and others (1953) is 
shown in figure 5. The agreement with the experimental curve, probably with- 
in the uncertainties of the two sets of data, suggests that a fairly close ap- 
proach to equilibrium has been attained in the experiments. 

In the second place, the experimental curve can be used to improve the 
thermochemical data, and thus to reduce the uncertainty in calculating the 
phase boundaries in other reactions involving albite, nepheline, and jadeite. 
The slope of the equilibrium curve is given by dP/dT = AS/AV, where AS 
is the difference in entropy between the products and the reactants, AV the 
difference in volume. Kelley and others (1953) calculate a slope of 20 
bars/deg; MacDonald (ms), in a critical review of the data, gives 20 + 3 
bars/deg. The required data are shown in table 3. These values are in such 
close agreement with the experimental value of 18 + 1 bars/deg that we can- 
not improve the values for AS or AV. It should be noted that, though in 
general dP/dT is a function of T(or of P) and the phase boundaries should 
show curvature, in the present study they are assumed to be straight. Neither 
the thermochemical nor the experimental data are of sufficient accuracy to 
show minor changes of slope of the equilibrium line, but both sets rule out 
changes greater than about 5 bars/deg for temperatures below 1200°C. 


3 
Entropies and Volumes, with Estimated Uncertainties 


Entropy Volume 

Phase Composition cal/deg mole cm*/mole 
a-quartz SiO. 10.0 + 0.1 22.64 + 0.03 
nepheline NaAlSiO, 29.7 + 0.3 54.2 +0.5 
jadeite 31.9 +05 60.7 0.6 
albite NaAlSi,0. 50.2 + 0.5 100.8 + 1.0 


The position of the equilibrium line depends upon the heat of reaction 
AH, the equation of the line being AG(T, P) = AH — TAS = 0, where G 
is the Gibbs potential. Table 4 shows the values of the heat of reaction for 
ab + ne = 2jd, according to the sources of the materials. The heat of this 
reaction is not measured directly but is obtained as a difference involving 
three relatively large numbers, the heats of solution of the three minerals. 
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Errors of 2-3 percent in the heats of solution could account for the whole 
range of values of AH. 

The calculations of Kelley and others (1953) depend upon the selection 
of AH = -5.5 kcal/mole, corresponding to the values for Varutrisk albite 
and Burmese jadeite; these are also the materials which they used in their 
studies of heat capacity. If we take the equilibrium pressure for the ab + ne 
= 2jd reaction, at 25°C, to be 1500 bars as found by extrapolation of the 
experimental line, we obtain AH = —3.6 + 0.8 kcal/mole, for this tempera- 
ture. Though lower than the average, and lower than the value adopted by 
Kelley and others, this is within the range given by Kracek and others (1951) 
and is probably not seriously inconsistent with their values. 


TABLE 4 
Heat of Reaction, ab + ne = 2jd, at 25°C, 1 bar 
(Kracek, Neuvonen, and Burley, 1951) 


H 
Albite Nepheline Jadeite kcal mole 

Varutrask Synthetic Burma -5.5 
Varutrask nis Japan, ground in agate mortar -6.7 
Varutrask . Japan, ground in mullite mortar -3.6 
Amelia Burma -7.2 
Amelia <j Japan, ground in agate mortar -8.4 
Amelia % Japan, ground in mullite mortar -5.3 


With this value, a revised calculation may be made for the reaction, 
albite = jadeite + quartz, which has been considered by Kracek and others 
(1951) and by Kelley and others (1953.) From the data of table 3, the slope 
of the equilibrium line is found to be 20 + 4 bars/deg. Let us suppose that 
the heat of the reaction, ab + ne = 2jd, is —3.6 + 0.8 kcal/mole and that 
some or all of the heats of solution for albite, nepheline, and jadeite are in 
error. If all of the error is attributed to the jadeite heat of solution, then the 
heat of reaction for ab = jd + qu becomes +1.0 kcal/mole. If all of the 
error is attributed to the heat of solution of albite, then the heat of reaction 
becomes +1.9 kcal/mole. If we take a mean value, +1.5 + 2 kcal/mole, 
the equilibrium curve for the breakdown of albite to jadeite plus quartz is 
given by P = 9000 (+5000) + 20 (+4) T, with P in bars, T in degrees 
centigrade. It is perhaps worth remarking that jadeite cannot be formed from 
nepheline plus quartz at relatively low pressures, as seems to be suggested by 
Kelley et al. In a mixture having the composition of jadeite, nepheline and 
quartz first react to form albite, according to ne + 2 qu = ab (Reaction C 
of Kracek et al.) at all points in the ab + ne field of figure 5. The albite will 
then react with the remaining nepheline to form jadeite only in the jadeite 
field. The reaction of nepheline and quartz to form albite was observed a 
number of times in the course of these experiments, but in only one case (no. 
61) did the reaction proceed far enough to produce jadeite. It is probably 
more advantageous to begin with albite and nepheline, to avoid the additional 
reaction. 
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_Fig. 7. Theoretical temperature-pressure relations in the plane of (dry) albite com- 
position. 

The P-T relations in the plane of albite composition are shown in figure 
7. At pressures somewhat above those of the albite breakdown, the quartz 
transforms to a high-density polymorph, “coesite” (Coes, 1953; MacDonald, 
unpublished). The slope of the albite melting curve has been computed from 
the change of entropy on melting (Kracek and Neuvonen, 1952) and the 
volume change (Goranson, 1938) as 40 bars/deg or 25°/kilobar; the un- 
certainty may be as large as + 40 percent. 

It is interesting to attempt the construction of temperature-composition 
planes at various pressures for the system NaAISiO,—NaAlSi,0, on the basis 
of the two P-T diagrams (figs. 5 and 7) of constant composition. Figure 8a 
shows the relations at one atmosphere as determined by Greig and Barth 
(1938.) Possible forms of the temperature-composition relations at 15,000 
bars and at 30,000 bars are shown in figure 8b and 8c. The liquidus relations 
are largely hypothetical but are somewhat controlled by the preliminary results 
of melting of material of jadeitic composition at 25,000 bars, and by the 
estimated increase of melting temperature of albite already cited. We have 
no experimental results as yet on the behavior of NaAlSiO, at high pressures; 
it seems likely, however, that nepheline will either break down to two or more 
phases of high density or transform to a high-density polymorph, indicated 
as nepheline - II in figure 8b and 8c. The possibility of such a form, having 
an olivine structure, is suggested by the relations in the system CaMgSi.0,— 
NaAlSiO,—SiO., for which a preliminary diagram by Schairer is given by 
Yoder (1950b, p. 318).° 


* Since the construction of figure 8, evidence for a transformation of nepheline at high 
pressures has been found by Mr. F. R. Boyd, Jr. at the Geophysical Laboratory. 
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Fig. 8. The binary system, NaAlSiO,—NaAlISi,Os, at several pressures. The diagram 
at 1 bar after Greig and Barth (1938). Liquidus relations largely hypothetical. 


Certain features of the postulated liquidus relations require comment. 
The eutectic is shown as moving toward the NaAlSiO, end of the system; 
this is suggested by the rapid increase of melting point of albite, but the 
amount of displacement is unknown. It is quite possible that an albite plus 
liquid field should be indicated in figure 5, and that the eutectic, at 15,000 
bars, has passed to the nepheline side of the jadeite composition. The ab + 
2jd equilibrium curve must reach a melting curve somewhere in the 


ne 


region indicated. We have supposed that beyond about 25,000 bars, jadeite 
melts congruently; this again requires verification. 

The significant feature illustrated in figure 8 is the reduction of the 
fields of stability of albite and nepheline. At 1 bar these minerals are stable 
from ordinary temperatures up to about 1000°C. At 30,000 bars the assem- 
blage albite plus nepheline is nowhere stable; the minerals nepheline and 
albite are stable only in the presence of jadite or liquid, except possibly for 
compositions near the end points at temperatures close to melting. 

Geologic considerations.—In applying these results to the Earth, a major 
restriction must be kept in mind. Only the system NaAlSiO,—NaAlSi,O, has 
been studied, and very few rocks can be adequately described in terms of 
these two components. The inferred albite breakdown, which would not be 
affected by the presence of excess silica, is applicable to a larger number of 
visible rock types. We probably do obtain, however, maximum pressures and 
roughly maximum depths at which various types of rock may have originated. 
When these observations are compared with estimates of internal tempera- 
tures, certain speculations as to the nature of the Mohoroviéié discontinuity 
begin to assume a new plausibility. 

Only a few natural occurrences of relatively pure jadeite have been 
described; besides those reported by Yoder (1950b), there are several new 
localities (Yoder and Chesterman, 1951; Wolfe, 1955). One of the classic 
localities in upper Burma has jadeite in association with albite, or nepheline, 
or both (Yoder, 1950b, p. 237). Our results indicate temperatures of crystal- 
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lization of 400-700°C if the pressure were 15,000 bars (with lower tempera- 
tures for lower pressures) in order to explain the jadeite-albite assemblage. 
Similiarly, high pressures and relatively low temperatures would be required 
for the albite-jadeite-quartz assemblage found in Japan (Yoder, 1950b, 
p. 240). 

The nepheline-albite assemblage is common in nepheline syenites. Both 
minerals, as found in nature, are in the form of complex solid solutions. The 
nephelines of plutonic rocks contain KAISiO,, SiO., and CaAl,Si,Og in addi- 
tion to NaAlSiO, (Miyashiro, 1951.) The feldspar may contain appreciable 
amounts of KAISi,O, and CaAl.Si.Og in addition to NaAlSi,O,. In order to 
discuss the stability of complex solid solutions at high pressures, we require 
data on excess volumes and excess enthalpies. MacDonald (ms) shows that 
an alkali feldspar containing both NaAlSi,O, and KAISi,O, will break down 
to jadeite plus orthoclase and quartz at lower pressures than those at which 
pure albite will break down. A similar effect may be anticipated in the case 
of other components than NaAlSiO, in nepheline, The experimental curve of 
figure 5 would then fix an upper limit for the pressure at which nepheline 
syenites could have formed, and with the liquidus would define a pressure- 
temperature region for formation of nepheline syenites. In nature this region 
would become more restricted, since the liquidus would be lowered by other 
components, especially water, and the equilibrium curve would be shifted 
toward lower pressures. 

In the same way, the line for the reaction, ab = jd + qu, fixes maxi- 
mum pressures (with due regard for the large uncertainty) for the formation 
of alkali feldspar or alkali feldspar-quartz assemblages. Such abundant rocks 
as granites and granodiorites probably crystallized at pressures well below 
those for this reaction. 

Eclogites present a more complicated problem. The pyroxene “omphacite” 
in eclogite is generally a complex solution of such components as Ca(MgFe) 
Si,0,, NaFeSi,O,, and NaAlSi,O,. Since diopside, hedenbergite, and acmite 
are stable at all temperatures up to the melting points at 1 bar, it is probable 
that the pressure needed to form omphacite is lower than the pressure needed 
to form jadeite. The pressures required to form an omphacite pyroxene in an 
eclogite undersaturated with SiO, would be lower, perhaps much lower, than 
those for the reaction ab + ne = 2jd; the omphacites in quartz-bearing 
eclogites would require pressures no greater than those for the reaction ab = 
jd + qu. It seems likely that, above some moderate pressure, there will be 
continuous solid solution between diopside and jadeite. 

The considerations reviewed above bear on the general problem of the 
stability within the Earth of the feldspars and feldspathoids. While it seems 
clear that these minerals must become unstable at depths of the order of some 
tens of kilometers, the discussion is complicated by the circumstance that 
we must consider multicomponent systems and by the uncertainty of the 
temperature-pressure relations in and below the crust. Several estimates of 
temperature versus depth (Poldervaart, 1955, p. 113) are shown in figure 9, 
together with the equilibrium curves which we have been considering. The 
line, P = 1000 + 18.5T for the ab + ne = 2jd equilibrium, has been trans- 
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Fig. 9. Equilibrium curves for several reactions versus depth, and estimated tem- 
perature-depth curves. A, B, and C represent different estimates for average continental 
conditions; D is an estimate for temperatures below average ocean of 4km depth, The 
quartz inversion after Yoder (1950a), the quartz-coesite curve after MacDonald (1956). 


formed to a temperature-depth curve, on the usual heavy-liquid assumption; 
in the typical continental section, the density has been taken as 2.8 to a depth 
of 33 km, 3.3 for greater depths. The line for the reaction, ab = jd + qu 
has also been drawn, as well as the quartz-coesite line, after MacDonald (1956) 
and the a-8-quartz curve, after Yoder (1950a). The curve of lowest tempera- 
ture (A) corresponds to an extreme concentration of heat sources within the 
crust; it intersects the jadeite curve at about 12 km. Curve B, a more reason- 
able distribution, is based on the assumption that one-half of the surface heat 
flow originates below the crust; this intersects the jadeite curve at about 
40 km. Curve C of the reference, intended to give an upper limit for a crust 
of normal thickness, might intersect the jadeite curve at some depth greater 
than 50 km, depending upon further assumptions about conductivity at depth 
and distribution of heat sources. 

Estimates for a typical oceanic section are shown in the small corner 
insert; the temperatures correspond to Case D of the reference, and the jadeite 
equilibrium line has been transformed to correspond to a structure composed 
of 4 km of water, 6 km of rock of density 2.8, with a density of 3.3 there- 
after. With these temperatures, points above a depth of about 37 km would 
lie in the albite-nepheline field, and the discontinuity at about 10 km below 
sea level could not be accounted for by the simple jadeite transformation. 
With allowance for additional components, however, the transformation will 
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probably shift to lower pressures, and none of these numbers should be 
regarded as finally disposing of the matter. 

Whether or not it proves possible to account for the Mohoroviéié dis- 
continuity in these terms, it seems safe to conclude that quartz, the feldspars, 
and the feldspathoids would transform to high-density forms at relatively 
moderate depths in the Earth. From the known properties of jadeite and 
garnets, we may conclude that these high-density forms would have elastic 
properties close to those of olivines and orthopyroxenes. The implications for 
geochemistry and petrology, and for our understanding of the structure and 
evolution of the Earth, are far-reaching. 


REFERENCES 
Adams, L. H., 1953, A note on the stability of jadeite: Am. Jour. Sct, v. 251, p. 299-308. 


Adams, L. H., and Gibson, R. E., 1929, The elastic properties of certain basic rocks and 
of their constituent minerals: Nat. Acad, Sci. Proc., v. 15, p, 713-724. 


Birch, Francis, 1939, Thermoelectric measurement of high temperatures in pressure ap- 
paratus: Rev. Sci. Instr., v. 10, p. 137-140. 


Bridgman, P. W., 1931, The physics of high pressure: London, G. Bell & Sons. 


1935, Measurement of certain electrical resistances, compressibilities and 
thermal expansions to 20,000 kg/cm*: Am. Acad, Arts Sci. Proc., v. 70, p, 71-101. 


1938, The resistance of nineteen metals to 30,000 kg/cm*: Am. Acad, Arts 
Sci. Proc., v. 72, p. 157-202. 


, 1940, The measurement of hydrostatic pressure to 30,000 kg/cm*: Am. Acad. 
Arts Sci. Proc., v, 74, p. 1-10. 


, 1948, The compression of 39 substances to 100,000 kg/cm*: Am. Acad, Arts 
Sci. Proc., v. 76, p. 55-87. 
Brothers, R. N., 1954, Glaucophane schists from the North Berkeley Hills, California: 
Am. Jour. Sct., v. 252, p. 614-626. 
Coes, Loring, Jr., 1953, A new dense crystalline silica: Science, v. 118, p, 131-132. 
Coes, Loring, Jr., 1955, High pressure minerals: Am, Ceramic Soc. Jour., v. 38, p. 298. 
Eskola, Pentti, Barth, T. F. W., and Correns, C. W., 1939, Die Enstehung der Gesteine: 
Berlin, Springer Verlag. 
Fermor, L. L., 1913a, Preliminary note on the origin of meteorites: Asiatic Soc. Bengal 
Jour. and Proc., v. 8, p. 315-324. 
—, 1913b, Preliminary note on Garnet as a geological barometer, and on an 
infra-plutonic zone in the earth's crust: Records, Geol. Survey of India, v. 42, p. 41-47. 


—, 1914, The relationship of isostasy, earthquakes, and vulcanicity to the earth's 
infra-plutonic shell: Geol. Mag., v. 1, p, 65-67. 


-, 1938, Garnets and their role in nature,: Indian Assoc. Cultivation Sci, Spe- 
cial Pub. 6. 


Goranson, R. W., 1938, Silicate-water systems: phase equilibria in the NaAlSisO;—H,O 
and KAISis0;—H.0 systems at high temperatures and pressures: Am. Jour, Sct., 5th 
ser., v. 35-A, p. 71-91. 


Greig, J. W., and Barth, T. F. W., 1938, The system NazO.Al.03.2Si0. (nephelite, car- 
negieite )—NasO.Al,03.6SiO. (albite): Am. Jour. Sct., 5th ser., v. 35-A, p, 93-112. 

Griggs, D. T., Kennedy, G. C., and Fyfe, W. S., 1955a, A simple apparatus for high pres- 
sure and temperature (abs.): Am, Geophys, Union Trans., v. 36, p. 511. 


Griggs, D. T., Fyfe, W. S., and Kennedy, G. C., 1955b, Jadeite, analcite, and nepheline- 
albite equilibrium (abs): Geol. Soc. America Bull., v. 66, p. 1569. 


Determination of Jadeite Stability Relations to 25,000 Bars 137 


Holmes, Arthur, 1945, Principles of physical geology: New York, The Ronald Press Co. 
Joly, John, 1930, The surface history of the earth: Oxford, Clarendon Press, 
Kelley, K. K., Todd, S. S., Orr, R. L., King, E. G., and Bonnickson, K. R., 1953, Thermo- 


dynamic properties of sodium-aluminum and potassium-aluminum silicates: U. S. Bur. 
Mines Rept. Inv, 4955. 


Kracek, F, C., and Neuvonen, K. J., 1952, Thermochemistry of plagioclase and alkali 
feldspars: Am. Jour. Sct., Bowen v., p. 293-318. 

Kracek, F. C., Neuvonen, K, J., and Burley, Gordon, 1951, Thermochemistry of mineral 
substances. I, A thermodynamic study of the stability of jadeite: Washington Acad. 
Sci. Jour., v. 41, p. 373-383. 

MacDonald, G. J. F., (ms.), 1954, A critical review of geologically important thermo- 
chemical data: Ph.D, thesis, Harvard University. 


, 1956, Quartz-coesite stability relations at high temperatures and pressures: 
Am. Jour, Sct., v. 254, p. 713-721. 


MacKenzie, W. S., 1954, Annual Rept. Director, Geophys. Lab., 1953-1954, p. 119. 

Miyashiro, Akiho, 1951, Ranges of chemical composition in nepheline and their petro- 
genetic significance: Geochim. Cosmochim. Acta, v. 1, p. 278-283. 

Poldervaart, Arie, ed., 1955, Physics of the crust: Geol. Soc, America Special Paper 62. 


Robertson, E. C., Birch, Francis, and MacDonald, G. J. F., 1955, Fields of stability of 
jadeite, kyanite, and pyrope (abs.): Geol. Soc. America Bull., v. 66, p. 1608. 


Thompson, J. B., Jr., 1955, The thermodynamic basis for the mineral facies concept: Am. 
Jour. Sct., v. 253, p. 65-103. 


Wolfe, C. W., 1955, Crystallography of jadeite crystals from near Cloverdale, California: 
Am. Mineralogist, v. 40, p. 248-260. 


Yoder, H. S., Jr., 1950a, High-low quartz inversion up to 10,000 bars: Am, Geophys. 
Union Trans., v. 31, p. 827-835. 


—_—, 1950b, The jadeite problem: Am. Jour. Sct., v. 248, p. 225-248, 312-334, 
—————., 1954, Annual Rept. Director, Geophys. Lab., 1953-1954, p, 121. 


Yoder, H. S., and Chesterman, C. W., 1951, Jadeite of San Benito County, California: 
California Div. Mines Special Rept. 10-C. 


Yoder, H. S., and Weir, C. E., 1951, Change of free energy with pressure of the reaction 
nepheline + albite — 2 jadeite: Am. Jour. Sct., v. 249, p. 683-694. 


[AMERICAN JOURNAL oF SCIENCE, VoL, 255, Fepruary 1957, P. 138-141] 


A NEW STANDARD FOR 
ESTIMATING DRAINAGE BASIN SHAPE 


RICHARD J. CHORLEY, DONALD E. G. MALM 
and HENRY A. POGORZELSKI 


ABSTRACT. The inapplicability of strict measures of circularity to descriptions of 
drainage basin shape is recognized, and one loop of a lemniscate of the form p= lcosk @ 
s suggested as a new standard, This allows for the development of two useful measures: 
k, which expresses a relationship between maximum basin length and maximum width 
which may have regional or structural implications; and the Lemniscate Ratio, which is 
a measure of how closely the actual drainage basin shape approaches that of its ideal 
lemniscate counterpart. 


The science of geomorphometry, which treats of the geometry of land- 
scape, takes as the basic unit of study the erosional drainage basin. Hitherto 
the shape of this unit has been expressed by the use of the standard methods 
employed in sedimentary petrology to describe the degree of rounding of 
particles, such as by comparing the area of the projection of the unit. with 
the area of a circle having the same perimeter as that of the unit ‘(Miller, 
1953). This method has two very obvious disadvantages, when applied to the 
shapes of drainage basins: 

1. Few, if any, drainage basins approach circularity, and there seems 
to be no tendency in a dissected region of more or less uniform structure for 
the steady state basin shape to be in any way related to the circle. Circularity, 
however, would seem to be an ideal end form for particles under many con- 
ditions of sedimentary transport, with which particle shapes can legitimately 
be compared. 

2. A numerical circularity ratio of the above sort gives very little visual 
impression of the actual shape of a drainage basin. 

It may be supposed that an “ideal” shape does exist for a drainage basin 
developed upon homogeneous bed rock, and this has been recognized by the 
common use of such terms as “pear-shaped” and “tear-shaped.” This ideal 
shape is capable of more precise definition as it has many properties of one 


petal or loop of a lemniscate (fig. 1), expressed by the equation: 
p 1 cos k 6 


Fig. 1. The lemniscate. 
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Where: lies between —7/2k and 7/2k. 

| is the longest diameter of the loop and in this instance equals 
the value of p when @ equals zero. It is analogous to the length 
of the drainage basin. 

k is a constant, greater than or equal to unity, upon which the 
relationship between the maximum length (1) and the maxi- 
mum width of the loop depends, — its control of the 
angle a. 

The value of k is interesting because its variations drastically control the 


shape of the loop (fig. 2). It should be noted that when k equals unity the 
loop is circular. 


Fig. 2. The relation between k and the lemniscate shape. 


The area (A) of the loop can be obtained in the following manner: 


/2k 
A= (l cos k 6)*. dé 


/2k 
4 2k —2/2k 
4k 
k = (1) 


As a first step in accurately describing drainage basin shape it is neces- 
sary to determine the lemniscate shape which the basin most. nearly ap- 
proaches. This can be accomplished by sttbstituting the measured values of | 
and A in equation 1 to obtain a value for k. 

Although this value, of itself, gives a good general impression of shape, 
it is profitable to further calculate how nearly the actual drainage basin shape 
approaches that of its ideal lemniscate counterpart. This can be accomplished 
by comparing the measured perimeter of the drainage basin (Pm) with the 
value of the perimeter of the ideal lemniscate (P) having the calculated value 


of k. 
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f (2 cos? k + Fk? sin?k 0). d 8 
e 2k 


on /2k 
21 V (cos? k 6 + sin? k 0). d 6 
0 


Let ¢6 = k @ such that when 6 = 0, ¢ = 0 and when 6 = 2/2k, ¢ = 2/2. 


V (cos? & + k* sin® 6).d¢ 
0 


Now let = 2/2-¥ 


a ¢7/2 
P= V (cos? (7/2 + sin? 
21 
=> V (sin? + k* cos? ¥).d¥ 
0 


°T 2 2 
= Vv (1 (k sin? ). dv 
e 0 k 


(k? — 1) 
Now if weletK = A aa perimeter (P) of the ideal lemniscate is 


k 


given by the equation: 


p=af V (1-—K?* sin? (2) 
0 


This is a complete elliptic integral of the second kind (usually denoted 
by E(K, +/2) or simply E(K)) and may be solved by the use of tables (Jahnke 
and Emde, 1945, p. 70-72), as follows: 

1. From the known value of k calculate K. 

From the tables find the value of arcsin K (i.e. sin~* K). 

From this, in the tables, find the value of E(K). 

Since: P = 21 E(K) (3) 
the value of the perimeter (P) can be easily calculated. 


It is now possible to calculate how closely the actual drainage basin shape 
approaches that of an ideal lemniscate loop by comparing the actual measured 


drainage basin perimeter with that of the ideal lemniscate, as a Lemniscate 
» 


I 
Ratio: p>: Obviously, as the value of this ratio approaches the upper limit- 


m 
ing value of unity, so the shape of the basin approaches that of the ideal 
lemniscate loop, i.e. that having the calculated value of k. 

Figure 3 illustrates how drainage basin shapes may be easily compared 
and emphasizes how true variations in shape are brought out by the use of 
the calculated k. Indeed, ranges of values of k may well be used to character- 
ize the regional tendency of drainage basin shape and thus be useful in dis- 
tinguishing one morphometric region from another. Another obvious applica- 
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A = 3.06 square miles A = 1.80 square miles 

1 = 2.00 miles 1 = 2.35 miles 

= 7.70 miles = 6.15 miles 
Circularity Ratio = Circularity Ratio = .598 
k = 1.026 k == 2.407 

K = .2238 K = .9094 

Arcsin K = 12° 56’ “t  Aresin K = 65° 25’ 
E(K) = 1.5509 E(K) = 1.1599 

P = 6.2036 P = 5.4515 

Lemniscate Ratio = .806 Lemniscate Ratio = .886 


Fig. 3. Calculation of circularity ratio, k, and lemniscate ratio for two drainage basins. 


tion lies in the use of k to express quantitatively structural control over basin 
shape, as, for example, in the effects of varying angles of dip on the shapes 
of drainage basins developed on a cuesta dip slope. The Lemniscate Ratio, 
too, has the possible merit of being used as an index of adjustment of drain- 
age basin shape to the ideal or “steady state” shape, as characterized by the 
calculated value of k. 
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X-RAY MEASUREMENT OF 
THE IRON-MAGNESIUM RATIO IN BIOTITES* 


JOHN ARTHUR GOWER 


ABSTRACT. The iron-magnesium ratio in biotite may be determined by measuring the 
intensity ratio of (004) to (005). This ratio is highly sensitive to iron-magnesium substi- 
tution and is virtually unaffected by other substitutions such as Na for K, F for OH, and 
Al for Si. The intensity method gives adequate precision and appears to be accurate to 
about 5 percent or less of the magnesium-iron content. 

Unit cell measurements are of limited value in the determination of compasition. ich 
rich biotites have a shorter co parameter than magnesium members, in spite of the fact 
that iron has a larger ionic radius than magnesium, Fluorine produces a marked decrease 
in co. Tetrahedrally coordinated Al may cause a slight increase. In metamorphic biotites 
co is larger than in igneous ones. The bo parameter shows little or no change with iron 
content but is smaller in fluorine-rich specimens, The a. parameter seems to show about 


the same variation as b. but the reflections in some specimens are too weak to be identified 
with certainty by powder methods. 


A new form of analysis plot has been used which shows more extensive solid solution 
with muscovite than has been previously supposed. 
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INTRODUCTION AND STATEMENT OF THE PROBLEM 


For many years the composition and abundance of feldspars have belie 
used to classify igneous rocks. The feldspars are useful minerals for this pur- 
pose because of the ease with which variations in composition may be deter- 
mined. However, feldspars tell only half the story of a rock. They say nothing 
of the ferromagnesian content and of the variation in relative amount of iron 
and magnesium. 

The amounts of the various ferromagnesian minerals in a rock may be 
easily determined petrographically. The relative amount of iron and mag- 
nesium in any of the ferromagnesians is not so easily determined, although in 
olivine, pyroxene, and to a certain extent amphibole, a rough estimate of the 
iron content may be made on the basis of optical properties. X-rays have also 
been useful in certain instances (Hess, 1952), where the iron content has been 
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related to the size of the unit cell. The commonest ferromagnesian in salic 
rocks and the ferromagnesian most widely distributed is biotite, and it would 
be desirable because of this to have some simple method of determining the 
Fe/Mg ratio in biotites. 

It has further been shown (Nockolds, 1947) that the ratio of magnesium 
to iron in igneous rock micas reflects very closely the magnesium-iron ratio 
of the rock as a whole. This is understandable because the known properties 
and compositions of phlogopites and biotites indicate a complete solid-solution 
series between magnesium-rich and iron-rich members (Winchell, 1935). 

To determine the ratio of magnesium to iron, various techniques have 
been described in the literature. Winchell and Winchell (1951) give curves 
showing changes in optical properties with change of composition. These can 
be used to give a rough approximation of composition, but in some cases the 
data are not reliable, especially for replacement in other than the Mg-Fe sites. 

The relation of color to composition has been studied by Hall (1941a). 
The triangular plot of FeO, MgO, and TiO, with respect to color shows that 
although color is most certainly controlled by iron and titanium content it is 
of limited value as a diagnostic feature. Blue-green biotites are low in titanium 
and rich in iron. Red or red-brown biotites may be rich in titanium but may 
have any Mg-Fe ratio. The plot might be more useful if the Fe/Mg content 
were known precisely. The color would then perhaps give an indication of 
titanium content. A further complication, not evaluated by Hall, is the effect 
of manganese, which presumably also produces a reddish-brown color. Ac- 
cording to Hall biotites do not in general contain sufficient manganese to have 
any effect upon the color, Geological occurrence would no doubt be of as- 
sistance in predicting manganese content. The effect of change in oxidation 
state of the iron is not known because of the lack of ferric-iron biotites, but 
for the analyzed specimens used by Hall it apparently has little effect on color. 

The relation between refractive index and composition has also been 
studied by Hall (1941b). Hall’s plot of refractive index (y) versus total iron 
shows that although iron causes an increase in refractive index, other elements 
affect it so much that correlation between iron content and refractive index 
is not possible. Titanium, manganese, and both ferric and ferrous iron raise 
the refractive index. Fluorine lowers the refractive index. 

Hall questions Winchell’s (1935) diagram relating chemical composition 
to optic properties because titaniferous biotites have much higher indices than 
even siderophyllite. The effect of 1% TiO. on refractive index is calculated by 
Hall to be + .0046, somewhat lower than the amount (.01) given by Kunitz 
(1936). 

The conclusion reached by Hall is that no sure information can be ob- 
tained as to the chemical composition of a biotite from its wefractive index 
alone, since a biotite high in iron and low in titanium may have the same in- 
dex as one low in iron and high in titanium. 

Heinrich (1946, 1953) states that the effect of ferric iron in raising re- 
fractive indices in biotite is approximately twice that of ferrous iron and 
about the same as that of titanium. Observations on synthetic fluor-phlogopite 
prepared by the Electro-technical Laboratory at Norris, Tennessee show that 
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replacement of all OH by F (about 9% by weight) lowers y by 0.22. This 
value is about half of the change in y produced by complete substitution of 


iron for magnesium. 


Taste | 
Structure Amplitudes (F) for Basal Reflections 


No. z 001 002 003 004 005 
B23 3.0 Me 0 300 203 275 25.7 236 
ry 12 40.9 31 -272 -$5,.2 -245 
28-74 365 205 266 -280 
6.0 O 34 -29.9 20.5 39.9 -210 -85 
1.0 K 5 -174 165 -149 182 -120 
Total 16.2 $.1 45.8 8.8 -49.4 
Feate % 8.0 10 4.4 -245 
Fos 5.2 25 206 N10 259 
AF 28 1.5 2.2 6.6 1.4 
B13 2.7 Mg 
tt 32.4 318 29 278 256 
0.03 li 
1.5 OH 
0.5 F 12 410.9 3.0 26.2 -33.4 -22.7 
4.0 O 
2.97 Si 
0.95 Al 23-75 -36.1 204 265 -281 
0.08 P 
6.0 O 34 9.3 20.5 40.5 21.0 —- 8.5 
5-16.05 15.25 -137 123 -121 
Total 19.6 ~ 6.6 58.9 12.2 —45.8 
% 90 - 3.0 26.9 5.6  -20.9 
!% 6.26 345 219 124 237 
AF 2.7 0.45 5.0 68 2.8 
B19 0.9 Mg 
1.15 Fe 
0 494 47.7 447 408 37,1 
0.3 Ti 
0.3 Al 
18 OH 
0.5 F 12 398 3.0 -263 -334 -228 
3.7 O 
or 28-74 -360 25 264 -276 
6.0 O 34 -293 -205 404 -210 -85 
0.74 K 
0.06 Na 50 -151 143 -129 17 -105 
0.09 Ca 
Total 37.4 85 664 245 -323 
Fresic % 14.2 3.2 25.2 93 -12.3 
Fos |% 18.7 36 192 130 169 
AF 15 0.4 6.0 3.7 4.6 


33.9 


- 44 


-13.7 
22.1 
9.4 
47.3 
17.9 


14.3 
3.6 


31.3 


005 007 
215 19.7 
-48 117 
13.9 295 
226 -14.4 
108 - 97 
36.2 368 
18.0 18.3 
19.3 15.5 
1.3 28 
23.5 213 
-138 293 
22.0 -14.0 
10.0 -89 
37.4 382 
7.1 175 
16.95 15.25 
0.15 2.25 
0.5 
29.1 
-14.0 
87 
47.2 
17.9 
14.3 
3.6 
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TABLE 1 (Continued) 


No. z 001 002 003 004 005 006 
B 7 0.68 Mg 
1.46 Fe+ + 
0.355 Fe 0 5387 52,06 48.56 4422 40.13 36.48 
0.225 Ti 
0.07 Al 
12 3980 3.02 -260 -32.7 -223 - 427 
28-74 -362 204 -280 --13,75 
6.0 34-292 -2045 40.5 -210 85° 22.05 
0.025 Ca 
0.130 Na 5S -1482 14.06 -12.7. 11.42 -1117 93 - 834 
0.7535  K 
Total 12.25 1249 69.76 28.44 -29.84 49.81 
Fesie % 49) 
Fis % 84 35 247 141:°°:170 199 
B9 2.05 Fe 
0.03 Mg 
0.35 Li 0 582 544 499 454 414 378 
0.07 Mn 
0.5 Al 
2.0 OH 
+ 12 410 31 -268 -35.1 -42 
28-74 -360 203 265 -278 -138 293 
6.0 0 34 295 -205 399 -210 -87 227 -144 
0.925 K 5 -169 162 -146 130 -118 106 -95 
Total 17 68.0 288 -313 526 511 


) 
53 23.2 98 -10.7 18.0 17.5 
Fons !% 8.5 25 26.7 14.5 55 25 183 
AF 7.0 2.8 3.5 4.7 4.8 2.5 2.2 


Optical spectrographic techniques may be used for semi-quantitative 
determinations of Mg and Fe, as well as for estimating relative amounts of 
Ti, Mn, Cs, Rb, Li, and many other trace elements. An important disadvantage 
lies in the fact that in studying rock biotites it is sometimes impossible to ob- 
tain a completely pure specimen. Admixed iron oxide, chlorite, amphibole, 
and intergrown crystals of zircon, rutile, sphene, and apatite would introduce 
large and unknown errors. These would exist not only in spectrographic 
analyses but also in standard chemical analyses. X-ray techniques obviate this 
last difficulty because small amounts of admixed minerals do not affect the 
structure of the mineral being studied. 


THE METHOD 
Development.—A common method of mineral identification using X-rays 
involves determination of lattice constants. For minerals in an isomerphous 
series such determinations are diagnostic if the lattice constants change 
markedly for the ion pair concerned but do not change for other composi- 
tional changes or for any other reasons. The changes in unit cell dimensions 
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are discussed later, but it is sufficient to say now that they are of limited use in 
determining the magnesium-iron content of biotites. 

Another possible method is to measure changes in intensity of X-ray re- 
flections with changes in magnesium-iron content. Brown (in Brindley, 1951) 
has calculated the variation in basal reflection intensities with variation in 
magnesium-iron content of biotites, and has suggested a method by which 
magnesium biotite may be distinguished from iron biotite using X-ray powder 
photographs. 

Because of varying degrees of preferred orientation in powder mounts 
and of different forms of layer stacking in the crystal structure it is desirable 
to use only basal reflections in the measurement of relative intensities, even 
though apparently one-layer monoclinic forms predominate (Heinrich and 
others, 1953). 

In order to determine the best pair of reflections to use for a precise 
determination of iron and magnesium content, the contribution of each atom 
to the amplitude of basal reflections has been calculated (table 1 and fig. 3). 
It was hoped that a situation exists where iron and magnesium make a large 
contribution to intensity in one reflection and practically none in another. 
Such is not the case. Although the contribution to intensity increases greatly 
with increase in iron content, the contributions to various OOl reflections 
change only slightly relative to one another. 

One significant feature is that the (005) reflections are opposite in phase 
to (004) reflections, whereas the Mg-Fe contribution to each of these is in the 
same phase. This means that as Fe is substituted for Mg the intensity of the 
(004) reflection increases and that of (005) decreases. Substitution for Si by 
Al, or for OH by F has practically no effect on intensities. Substitution of Na 
for K might interfere because of the difference in atomic number, and thus 
in scattering effect. Quite fortuitously the phases of the structure amplitude for 
K in (004) and (005) are opposite so that substitution of Na for K has a 
similar effect on both reflections and thus produces very little change in the 
ratio. 

Replacement of Mg or Fe by other cations is not extensive, although it 
may be enough to interfere with the curve. Manganese has almost the same 
scattering factor as Fe; Al has almost the same scattering factor as Mg. 
Lithium, having a very small scattering factor, would act in the same direc- 
tion as Mg but more strongly. Where these elements are greatly enriched, other 
minerals in the rock will frequently indicate their presence—e.g., spodumene, 
topaz, and muscovite. Where Al is abnormally low, olivines and pyroxenes 
may be associated with the biotite. 

Two types of working curve were prepared. The first (fig. 1) is a plot 
of the ratio of structure amplitudes (F’s) of (004) to (005) against atomic 
Fe 
Mg+Fe+Al 
points representing the analyzed specimens is straight. This is to be expected 
because F is a summation, one of whose terms contains the scattering factor 
(/). The second curve (fig. 2) shows the ratio of uncorrected intensities for 


proportion of iron [ 100 x ( )] . A line drawn through the 
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FeKa radiation against molecular percent of Fe, as above. This form is more 
readily used but because of a shift in glancing angles with compositions is 
less accurate than figure 1. 

Limitations—The use of the curves (figs. 1 and 2) is limited to biotites 
which contain principally magnesium and iron in the octahedral positions. 
The ratio actually measured depends on the relative amounts of Fe, Mn, Ti, 
Al (octahedral), Li, and Mg in the structure. However, the effects of Mn and 
Ti (f ~ f for Fe) tend to cancel those of Al and Li (f ~ f for Mg), and in- 


as much as high concentrations of these elements may frequently be detected 


‘Fig: Change ih Ratio ot Stracural Arhptitudes ot (006) to Q06)/with Change.’ 
thes 
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optically some adjustment could be made to the observed ratio to correct for 
them. 

A rare alkali biotite (Hess and Stevens, 1937) serves as a test of such a 
variety. In spite of the unusual composition of the mica the observed ratio, 
Fe/Mg+Fe %, falls very close to the curve (fig. 1). 

Reproducibility ——A high degree of precision is attainable if a great deal 
of care is exercised in preparation of the specimen. The source of greatest 
error lies in not having the surface of the powder perfectly flat and in the 
plane of the holder. An uneven surface produces irregularly shaped peaks 
which may give large errors in intensity ratios. 

The observed extent of these errors is indicated by vertical lines (figs. 1 
and 2) through the points representing the best values for the analyzed 
specimens. 

Sample preparation and measuring techniques——The method found to 
give the highest degree of precision is that described by McCreery (1949). 
The mica is first reduced to powder by filing the edges of cleavage pieces or, 
if the grains are too small to file, by grinding in acetone with 120 grade 
carborundum in an iron mortar. Impurities may be removed by heavy liquids 
and by magnetic separation. Although —400 mesh material gives the flattest 
surface coarser sizes (up to —270 mesh) are satisfactory especially if the mica 
is not well crystallized. A greater degree of preferred orientation is produced 
and basal reflections will thus have greater intensity. 

The powder is packed gently into a drilled 14” diameter hole in a thin 
section slide, backed by another slide. A third slide is placed on top and the 
whole turned over. The original backing slide is carefully lifted off so that the 
surface of the powder is not disturbed. The completed mount is then bound 
at one end with plastic tape or drafting tape. The latter was found more satis- 
factory both because of ease in labeling the specimen and ease in removal if 
the mount is to be reversed in the gonimeter. An internal standard may be 
mixed with the powder if desired. 

A North American Philips X-ray diffractometer and Brown recorder were 
used in all the work. Both copper and iron radiation were used, but iron was 
preferred because of improved patterns for iron-rich biotites. Preliminary 
traverses were made at 50 kv and 10 Ma at 1° per minute using 1° slits and 
a time constant of 2 seconds. This revealed any impurities, such as chlorite, 
amphiboles, or pyroxenes, and gave the approximate 26 angles for all reflec- 
tions. 

The (004) and (005) peaks were traversed at 14° per minute. Accuracy 
was improved by measuring a number of sample mounts a few times each 
rather than measuring the same mount many times. If the specimen is well 
crystallized 1° slits may be used, but if not, 4° slits are necessary to give peaks 
large enough to measure. 

Areas under peaks are measured with a planimeter, with which precision 
of 1 percent is obtainable by these repetitions of each measurement. 

Accuracy.—The accuracy of the method could not be fully investigated 
because of insufficient analyzed material. However, a spectrochemical method 
gives a rough check. A curve of line intensity ratio versus iron concentration 
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was plotted using the available analyzed specimens and the indicated iron 
content of several biotites was compared to that given by the X-ray method 
(fig. 7). Adequate correlation was obtained where spectrochemical determina- 
tions gave good reproducibility (about 1% spread in measured Fe/Mg+Fe 
ratio), but some difficulty was experienced in developing a technique for 
arcing biotites. X-ray determinations on specimens from localities where 
similar material had been analyzed also gave satisfactory agreement (speci- 
mens B 18, B 26, B 7). 


Fig. 3 
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UNIT CELL DETERMINATIONS 

Unit cell sizes were measured by using (004), (005), (060), and where 
possible (200). 

The powdered sample was mixed with an internal standard of finely 
powdered (-400 mesh) silicon and X-rayed as described above. Results for c 
sin B(doo:), 6, and a are given in table 2 and figure 4. 

For the chemically analyzed specimens and for those checked spectro- 
chemically a decrease in c sin B is noted with increase of iron content. At first 
sight this seems anomalous because the ionic radius of Fet+? (.74 A, Ahrens, 
1952) is larger than that of Mgt* (.67 A). However, similar changes are 
shown by chlorite (von Engelhardt, 1942). This phenomenon may have the 


| 
‘ 
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following explanation. Iron has more than twice the polarizing power of mag- 
nesium,’ and may so polarize the OH layer as to reduce the distance between 
apical oxygen atoms and thus partially collapse the structure. No great in- 
crease in b was found such as was found in chlorite by von Engelhardt. 


Fig.5 Variation in Fe,Mg, and 
Octahedral Al in Micas 


Al 9 


Using ze 
. + , polarizing power for Mg 335 and Fe = 720 (see Evans, 1948, p. 181, 


and Ahrens, 1953). 
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The effect of F was studied in a few analyzed phlogopites. Complete re- 
placement of OH by F reduces c sin 8 by about 0.1 A (see fig. 4) and a 
rough measure of F content may be made as indicated in figure 4. 

In iron-biotites where more aluminum commonly substitutes in both octa- 
hedral and tetrahedral layers, the effect of F is obscured—as for example in 
B 9, which is rich in both F and Al. Aluminum apparently increases c sin B, 
but lacking analyzed material positive tests could not be made. 

A test on metamorphic biotites from Vermont® indicate that they have a 
larger doo, spacing than do igneous. This may be caused by a high percentage 
of aluminum, a low confining pressure during crystal growth, a partial substi- 
tution of water for potassium, or some other cause. A granite with “cold” 
contacts from the same area contains a biotite with a larger-than-normal doo, 
spacing. Further work is required on this phase of the subject. 


TABLE 2 
co sin B, b,, and a, for Biotites 
Spec. no. c. sin B b. fo Fe/Fe Mg Al 
B 1 10.057 A 9.250 A 59 % 
2 10.032 9,223 33 
3 10.081 9.245 50 
7 10.046 9.250 71 
8 10.010 9,232 1 
9 10.024 0.239 89 
1] 10.073 9,242 5.23 55 
12 10.053 9,237 5.28 59 
13 10.076 9.235 5.23 9 
14 10.012 ~~ 92 
15 10.037 9.231 5.26 47 
16 10.103 9.254 30 
18 10.051 9.213 5.22 11 
19 10.051 9.243 5.28 61 
20 10.072 9.253 5.24 68 
21 10.072 9.250 5.24 54 
22 10.065 9.255 5.26 81 
23 9.985 9.179 5.22 0 
24 10.017 ~- 32 
25 10.087 9.247 5.28 76 
26 10.093 9.238 3 
27 10.042 9.231 55 
31 9.915 9.066 5.14 Muscovite 
34 10.066 9.216 51 
BeF 19 10.092 9.257 23 
Be 74 10.103 9.257 35 
G 43 10.095 9.251 29 
BeF 14 10.096 9.254 5.26 ? 53 
BH 9 3 10.095 9.260 48 


E 15 4 10.097 9.251 43 


? Lent by Professor G.J.F. MacDonald of the Department of Geology and Geophysics, 
Massachusetts Institute of Technology. 
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RELATION OF PHYSICAL PROPERTIES AND COMPOSITION TO ENVIRONMEN1 

Several studies have been made of the variation in composition of biotites 
with respect to variation in their parent rock type (Nockolds, 1947; Heinrich, 
1946, 1953). A few of the more significant conclusions are outlined below. 

In a study of Caledonian igheous rocks from Western Scotland, Nockolds 
(1947) found that in cale-alkali igneous rocks those biotites associated with 
highly aluminous minerals are rich in Al.O,, those associated with hornblende, 
pyroxene, or olivine are low in Al.O,, and those which occur alone are 
intermediate. He concludes that the Al/Mg+Fe ratio is determined by the 
paragenesis of the rock, but the Fe/Mg ratio depends on the degree of differ- 
entiation of the souree magma, or upon the amount of contamination by as- 
similation of country rock. In other words the Al/Mg+Fe ratio reflects the 
environment, whereas the Fe/Mg ratio reflects the composition of the rock 
as a whole. 

The explanation of this relationship may be as follows: Aluminum can 
enter into limited substitution with Mg and Fe in the octahedral positions of 
biotites. In aluminum-rich rocks, then, the biotites can take only part of the 
aluminum, and separate aluminum-rich minerals will be formed—such as 
muscovite, spodumene, and topaz. In aluminum-poor rocks the available 
aluminum will be used by biotite and the associated minerals will be non- 
aluminous, e.g., pyroxene, olivine, and to a certain extent amphibole. Substi- 
tution of Mg for Fe in biotites is complete, however, and the ratio of Mg to 
Fe in any biotite is probably close to that of the rock as a whole at the time 
of formation of the biotite. 

Heinrich (1946) studied the relation of the chemical variation of biotites 
to geologic occurrence. He plotted ratios of FeO+MnO, Fe.0,+TiO., and 


Fig 6 Variation in |ron Content of Biotites with Rock Type 
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MgO in biotites from eight rock types and found that the ratios fell into 
restricted fields. The ratio of FeO to Fe,O, does not vary greatly with rock 
type but the ratio of total Fe to Mg varies from close to zero to a very large 
figure equivalent to the sequence of Bowen’s fractional crystallization series. 

The explanation of this may lie in the relative bonding energies of mag- 
nesium to oxygen and iron to oxygen. Because of its greater ionization poten- 
tial magnesium tends to form a stronger bond with oxygen than does iron, so 
that magnesium silicates have in general a higher melting point than their iron 
analogues (Ahrens, 1953). 


CONCLUSIONS AND SUGGESTIONS FOR FUTURE RESEARCH 


The foregoing method of determining the extent of magnesium-iron sub- 
stitution in biotites by X-ray diffractometry makes use of a ratio of X-ray 
intensities that is very sensitive to changes in the amounts of magnesium and 
iron and relatively insensitive to other compositional changes. The chief dis- 
advantage of the method is that manganese, lithium, aluminum, and titanium, 
which may substitute in the octahedral sites, are not detected. 

Careful measurements of cell dimensions indicate that these are of little 
value in determining composition but do show systematic changes for which 
explanations can be suggested. 

The intensity method is rapid and reasonably accurate. It may find con- 
siderable use in studying changes in the magnesium-iron content of rocks in 
metamorphic zones, and in more precise classification of igneous rocks than is 
at present possible by petrographic methods. 

Further research, using similar methods of determining iron/magnesium 
ratios, is recommended for other minerals. In fact the method may be applied 
to any solid-solution series in which the substituting elements have a large 
difference in atomic number (and thus in atomic scattering factor). Some of 


Variation in’ Spactrateintensities with Composition. of | Jiotite 
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the more obvious series are: olivines, amphiboles, pyroxenes, carbonates, 
chlorites, montmorillonites, epidotes, garnets, spinels, and tourmalines. 


APPENDIX 


Chemical analyses were available for only five of the specimens studied. 
These are presented in table 4. Hundreds of biotite analyses are to be found in 
the literature. Compilations have been made by Doelter (1917), Hall (194la, 
1941b), Heinrich (1946), and Nockolds (1947). Many of them show poor 
agreement with ideal structural formulae and may have been made on impure 
material or by inaccurate methods. One hundred of the more reliable have 
been plotted on a triangular diagram (fig. 5) showing ratios of Fe, Mg, and 
octahedral Al. The diagram shows the biotite-phlogopite field, the muscovite- 
ferrimuscovite-picrophengite field, and a supposel “forbidden” field. However, 
some biotites apparently contain enough octahedrally coordinated aluminum 
to invade this field. It is noteworthy that all of these are associated with alumi- 
nous minerals such as topaz, spodumene, muscovite, and beryl. The source 
rocks for various biotites are indicated on the plot and occupy fairly definite, 


TABLE 4 
Chemical Analyses 


B 23 B 13 B 19* B7 B 9 

42.50% 42.19% 34.10% 33.07% 37.01% 
TiO. 0.61 4.90 3.84 0.02 
A].0 12.35 12.60 16.65 16.32 15.89 
Fe.0, 0.04 0.33 3.86 5.97 Trace 
FeO 2.91 17.75 22.46 30.16 
MnO 0.02 0.59 1.01 
McO 28.36 25.61 7.91 5.85 0.22 
CaO 0.14 0.06 1.08 0.26 0.10 
Na.O 0.06 0.31 0.45 0.87 0.58 
K.O 10.88 9.91 7.50 7.92 9.02 
Li.O 1.01 
0.19 
0.12 
F 9.30 2.66 2.04 3.88 
Cl 0.24 
H.0 + 2.84 3.48 3.87 1.92 
H.O 0.33 0.47 
P.O 0.13 
Ign. loss at 0.15 
640° C 

103.78 100.51 100.78 100.43 101.37 
Less O for F 3.92 1.12 0.89 — 1.68 

99.86 99.39 99.89 100.43 99.69 
* This specimen contains 8% chlorite (ny = 1.630, nz - nx = 0.00) for which adjust- 


ment is made. 
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although overlapping, fields. A greater tolerance for octahedral aluminum is 
indicated in iron-rich varieties. To show the iron-magnesium ratios directly 
the data are also presented in a more simplified plot (fig. 6) showing the 
range in iron-magnesium ratios for biotites of each major rock type. 

The desirability of determining the iron-magnesium ratio is clearly 
shown. Each rock type contains biotites of limited compositional range. It is 
important to know not only that the biotite composition falls in the range 
given but also what the composition is in terms of iron and magnesium. 
Mineral association and optics will indicate unusual enrichments of rare alka- 
lies, fluorine, and aluminum but will not give iron and magnesium content. 
The X-ray intensity method does this satisfactorily enough to be of some use 
in petrologic studies, 


List of Specimens 
Bl ___ Biotite from granite, Chester, Nova Scotia. 
B2___Biotite from gabbro, Tahawus, New York. 
B3__Biotite from Capilano quartz diorite, Vancouver, B.C. 
B7___—sCBiotite from Miask, Ural Mountains. 
B8  Phlogopite, Hadderspell Twp, Quebec. 
B9 Siderophyllite, Brooks Mountain, Alaska. 
Bll Biotite from alnoite, Storanset, Alno, Sweden. 
Biotite, unknown. 
B13 Phlogopite, west of Dent, Clearwater County, Idaho (U.S.G.S., A. Hietanen, personal 
communication ). 
B14 Biotite, associated with wohlerite, Brevig, Norway. 
B15 Biotite, Renfrew, Ontario. 
B16 Biotite, Asbestos mine southwest of Mt. Lincoln, Belchertown, Mass. 
B18 Phlogopite, Burgess, Ontario. 
B19 Biotite from “G 1” granite, Westerly, Rhode Island. 
B20 Biotite from granite, Brattleboro Quadrangle, Vermont. 
B21 Biotite from granite, West Nictaux, Nova Scotia. 
B22 Biotite from granite, Albany Cross, Nova Scotia. 


B23 Synthetic fluor-phlogopite, Electrotechnical Laboratory, Norris, Tennessee (R. A. 
Hatch, personal communication). 


B24 Biotite, Kings Mountain, North Carolina (ref, Hess and Stevens, 1937, p. 1045). 

B25 Biotite from Nordmarkite (?), Brome Mountain Quarry, Iron Hill Road, Quebec. 

B26 Phlogopite from dolomite, Simplon Tunnel (ref. Jakob, 1928, p. 219). 

B27 Biotite from essexite, Mt. Johnson, Quebec. 

B28 Synthetic phlogopite from Geophysical Laboratory, Washington (950°C, 15,000 
lb/in® water vapor pressure, 5 hrs.) (ref. Yoder and Eugster, 1954, p. 157). 

B31 Muscovite—labeled annite from Cape Ann granite, Rockport, Mass. 

B32 Biotite from syenite, Kingston Mills, Ontario. 

B33 Biotite from foot of Mt. Megantic, Quebec. 

B34 Phlogopite, Sebastopol, Ontario. 

Metamorphic Biotites from Vermont 

BeF 19 With garnet, staurolite (with chlorite alteration), quartz. 

Be 74 With kyanite, staurolite, chlorite. 

G43  Chlorite-biotite. 

BeF 14 Biotite porphyroblasts with chlorite rims, quartz, muscovite, and paragonite, 

BH 9 3 Garnet zone—with orthoclase. 

E 15 4 Garnet zone—no feldspar. 
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The Geology of the Bennington Area, Vermont; by Joun A. MACFADYEN, 
Jn. P. 72; 3 figs., 8 pls. Vermont Geological Survey Bulletin 7, Montpelier, 
Vermont, 1956 (Vermort Development Commission).—This bulletin is one 
issue of the new series whose attractive and readable format is becoming in- 
creasingly familiar to those concerned with the geology of Vermont. A notable 
feature of the text is the six-page section on the Taconic thrust problem; like 
others who have worked in the vicinity, the author rejects the hypothesis of 
an extensive klippe in eastern New York and western New England. 

The paper is a significant contribution to the descriptive geology of the 
south end of the Green Mountain anticlinorium and of the eastern part of the 
synclinorium to the west. The conclusions arrived at depend inherently upon 
data that must come from a region much larger than that of the Bennington 
area. Problems such as that of the relative importance of Taconic versus 
Acadian orogeny in the formation of the Green Mountain anticlinorium, or 
that of the Taconic thrust, are quite obviously not apparent within the confines 
of the geologic map of the Bennington area. The parable of the three blind 
men and the elephant seems quite applicable to this and similar situations, and 
in the reviewer's feeling points up the fact that understanding depends fully 
as much on judicious extensive investigation as upon more local studies that 
seem to supply “open and shut” information. 

The logic involved in the discussion of the age of the Grahamville forma- 
tion as reflecting the age of the Hoosac formation (p. 30) is to a considerable 
extent circular inasmuch as Thompson (1952, Geol. Soc. America Guidebook 
for Field Trips in New England, p. 18) based his inference as to the age of 
the Grahamville chiefly upon the apparent stratigraphic relationships of the 
Hoosac at its type locality between North Adams, Massachusetts and Heart- 
wellville, Vermont. WALLACE M. CADY 


Geologic Map of the State of Oregon, 1/500,000; by L. F. Hintze. Provo, 
Utah, 1956 ($20.00) —This, the first geologic map of Oregon, has been‘com- 
piled by L. F. Hintze, a consulting geologist, from U. S. Geological Survey 
open file material and other sources. The map represents Cenozoic geology 
faithfully, and it illustrates such features as the distribution of volcani¢ cones 
and high angle faults clearly and in commendable detail. The geology of areas 
of older rocks as in the Klamath and the Blue Mountains is more generalized, 
as metamorphic and sedimentary rocks are not differentiated. 

The map is an “ozalid” and must be hand-colored. The use of patterns 
for some rock units might have helped and might have been expected in justi- 
fication of the price. As the map does not show topography (except for 
hachured volcanic cones), the inclusion of more geographic information would 
have been useful. Although the map is welcome for the substantial amount of 
information set forth, we might reasonably expect a closer correlation between 
price and utility. LAWRENCE LUNDGREN, JR. 
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University Press, $7.50). 

Introduction to Microfossils; by D. J. Jones. New York, 1956 (Harper & Brothers, $6.50). 

The American Arbacia and Other Sea Urchins; by E. B, Harvey. Princeton, 1956 (Prince- 
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and Canada). Joint Hydrology-Sedimentation Bull, 7, Washington, 1955 (Superinten- 
dent of Documents, $1.25). 

North Dakota Geological Survey Rept. of Investigation: 22, North Dakota’s Nesson Anti- 
cline; by W. M. Laird and C. B. Folsom, Jr. 23, The General Geology of Uranium in 
Southwestern North Dakota; by J. R. Bergstrom, Giand Forks, 1956. 

Illinois Geological Survey: Rept. of Investigations: 194, Groundwater Geology of Lee and 
Whiteside Counties, Illinois; by J. W. Foster. 195, Preliminary Report on Portland 
Cement Materials in Illinois; by J. E. Lamar and others, Circular 216, Germanium 
in Fly Ash and Its Spectrographical Determination; by J. S. Machin and Juanita 
Witters. Urbana, 1956. 

Vermont Geological Survey Bulletin 8: The Geology of the Lyndonville Area, Vermont; 
by J. G. Dennis, Montpelier, 1956. 

1955 Vacuum Symposium Transactions, Boston, 1956 (Committee on Vacuum Techniques, 
Inc., Box 1282, $10.00). 


Geology and Ourselves; by F. H, Edmunds. New York, 1956 (Philosophical Library, 
$10.00). 


Microscopic Petrography; by E. Wm. Heinrich. New York, 1956 (McGraw-Hill Book 
Company, $6.50). 

Report of the Special Committee on the Federal Loyalty-Security Program of the Associa- 
tion of the Bar of the City of New York. New York, 1956 (42 West 44th Street, 
$5.00). 

Journal of the Franklin Institute Monograph 2: Earth Satellites as Research Vehicles. 
Philadelphia, 1956 ($2.50). 

The Generation of Electricity by Wind Power; by E. W. Golding. New York, 1956 (Philo- 
sophical Library, $12.00). 

The Petroleum Refinery Engineer's Handbook; by J. F. Strachan. New York, 1956 (Philo- 
sophical Library, $15.00). 

Action des Rayonnements de Grande Energie sur les Solides; by Y. Cauchois and others. 
Paris, 1956 (Gauthier-Villars, 1800 francs, paper cover; 2100 francs, cloth bound). 

Geology and Ground-water Resources of Rawlins County, Kansas; by K. L. Walters. 
Kansas Geol. Survey Bull, 117. Lawrence, 1956. 

Ground Water in Wisconsin; by W. J. Drescher. Wisconsin Geol. Survey Information 
Cire. 3, Madison, 1956. 

Bibliography of Seismology, July-December 1954; by W. E. T. Smith. Ottawa, 1956 
(Dominion Observatory, $ .25). 

Valency and Molecular Structure; by E. Cartmell and G. W. A. Fowles. New York, 
1956 (Academic Press, $5.80). 

Classification of the Pennsylvanian Rocks of Illinois as of 1956; by H. R. Wanless. Illinois 
Geol. Survey Circ, 217. Urbana, 1956. 

The Trickster; A Study in American Indian Mythology; by Paul Radin. New York, 1956 
(Philosophical Library, $6.00). 

Abacs or Nomograms; by A. Giet. New York, 1956 (Philosophical Library, $12.00). 

The History of the British Flora by H. Godwin. New York, 1956 (Cambridge University 
Press, $16.50). 

Physics, Psychology, and Medicine; by J. H. Woodger. New York, 1956 (Cambridge 
University Press, $1.75). 

Introduction to Scientific German, 2d ed.; by J. H. Wild. New York, 1956 (Oxford Uni- 
versity Press, $2.75). 


California and the Southwest; C. M. Zierer, editor, New York, 1956 (John Wiley & Sons, 
$11.25). 
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Measurements of Mind and Matter; by G. W. Scott Blair. New York, 1956 (Philosophical 
Library, $4.50). 

Minerals and the Microscope; by H. G. Smith. Fourth ed., revised; by M. K. Wells (New 
York, 1956 (Macmillan Company, $3.00). 

Dead Towns and Living Men; by Sir Leonard Woolley. New York, 1956 (Philosophical 
Library, $6.00). 

Elementary Crystallography; by M. J, Buerger. New York, 1956 (John Wiley & Sons, 
$8.75). 

Man’s Role in Changing the Face of the Earth; W. L. Thomas, Jr., editor. Chicago, 1956 
(University of Chicago Press, $12.50). 

The Earth Beneath Us; by H. H, Swinnerton. Boston, 1956 (Little, Brown and Company, 
$5.00). 

Ceramics for the Archaeologist; by Anna O, Shepard, Carnegie Inst. Washington Pub. 
609. Washington, 1956 ($6.75, paper bound; $7.75, cloth bound). 

The Potentials about a Point Electrode and Apparent Resistivity Curves for a Two-, 
Three- and Four-layer Earth; by H. M. Mooney and W. W. Wetzel. (With accompany- 
ing set of master resistivity curves in paper cover.) Minneapolis, 1956 (University of 
Minnesota Press, text only, $4.90; set of curves only, $15.00; complete, $18.00). 

Stratigraphie Griechenlands; by C, Renz, Athens, 1955 (Institute for Geology and Sub- 
surface Research, $12.00 plus mail charges). 

Wisconsin Geological and Natural History Survey Bull. 80: Soil Survey of Grant County, 
Wisconsin, an Introductory Report; by F, D, Hole. Madison, 1956. 

Kansas Geological Survey Bull. 119, Pt. 4: The Mineral Industry in Kansas, 1950-1954; 
by W. H. Schoewe. Lawrence, 1956. 

Georgia Dept. of Mines, Mining and Geology Circ. 2: Directory of Georgia Mineral Pro- 
ducers, 1956 (9th ed.). Atlanta, 1956. 

Illinois Geological Survey Rept. Investigations 196: Groundwater Geology of White County, 
Illinois; by W. A, Pryor, Urbana, 1956. 

Une Tentative d’Interprétation Causale et non Linéaire de la Mécanique Ondulatoire; by 
Louis de Broglie. Paris, 1956 (Gauthier-Villars, 3500 francs). 

Bird and Butterfly Mysteries; by Bernard Acworth. New York, 1956 (Philosophical Li- 
brary, $7.50). 

Geological Survey of Great Britain Bull, 11. London, 1956 (Her Majesty's Stationery 
Office, 17s. 6d.). 

Handbook of Scientific and Technical Awards in the United States and Canada, 1900- 
1952: Margaret A, Firth, editor. New York, 1956 (Special Libraries Association, 
$10.00). 

North Dakota Geological Survey: Rept. of Investigations: 24, Mississippian Stratigraphic 
Studies, Bottineau County; by S. B. Anderson and L. B. Nelson. 25, Geology of the 
Lower Pipestem Creek Area North Dakota; by R. J. Kresl. 26, The Geology of the 
Grassy Butte Area, McKenzie County, North Dakota; by E, G. Meldahl. Grand Forks, 
1956. 

Annual Report of the Geological Survey Department of Nigeria 1954-55, Lagos, 1956. 

Canada Dept. Mines and Technical Surveys Mem. Series 132: Interim Report Hardness of 
Major Canadian Water Supplies; by J. F. J. Thomas, Ottawa, 1956 ($ .25). 

Buffalo Society of Natura] Sciences Misc, Contributions: 2, Geology of Chautauqua 
County, New York; by LH. Tesmer. 3, Panorama of Summer Wild Flowers; by C. 
A. Zenkert, Buffalo, 1956. 

The Barber County Earthquake of January 6, 1956; by Louis F. Dellwig. State Geological 
Survey of Kansas, Bull, 119, Pt. 5, Lawrence, 1956. 

Geologia e Petrologia do Arquipélago de Fernando de Noronha; by Fernando F. M. de 
Almeida, Departamento Nacional da Producgéo Mineral Monografia XIII, Rio de 
Janeiro, 1955. 

University of Tasmania Department of Geology Publications: 40, Cyclic Sedimentation 
in the Permian System of Tasmania; by K. G. Brill, Jr. 41, The Geology of the Up- 
per Huon-Arve River Area; by R. J. Ford, 42, The Geology of the West Coast Range 
of Tasmania—Part I]—Structure and Ore Deposits; by J. Bradley. Hobart, 1956, 

Cellular Aspects of Basic Mechanisms in Radiobiology; H. M. Patt and E. L. Powers, 
editors, Nuclear Science Series Report No. 18. Washington, D, C., 1956. (National 
Academy of Sciences-National Research Council Pub, 450). 
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